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60/257624, "Intrinsic Phase Shifter as an Element of a Superconducting Phase 
Quantum Bit", filed December 22, 2000, from which priority is being claimed, 
herein incorporated by reference in its entirety. This application is further related 
15 to concurrently filed application Serial No. [Attorney Docket No. M-8915 US] 
entitled "Quantum Bit with a Multi-Terminal Junction and Loop with a Phase 
Shift" and application Serial No. [Attorney Docket No. M-8915-1C US] entitled 
"Quantum Bit with a Multi-Terminal Junction and Loop with a Phase Shift", both 
of which are herein incorporated by reference in their entirety. 
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BACKGROUND 


1. Field of the Invention 

This invention relates to quantum computing and, more specifically, to 
2 5 solid state quantum computing qubits with superconducting materials. 


2. Discussion of Related Art 

Research on what is now called quantum computing traces back to Richard 
Feynman, See, e.g., R. Feynman, Int. J. Theor. Phys., 21, 467-488 (1982). 
3 0 Feynman noted that quantum systems are inherently difficult to simulate with 
conventional computers but that observing the evolution of a quantum system 
could provide a much faster way to solve some computational problems. In 
particular, solving a theory for the behavior of a quantum system commonly 
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involves solving a differential equation related to the Hamiltonian of the quantum 
system. Observing the behavior of the quantum system provides information 
regarding the solutions to the equation. 

Further efforts in quantum computing were initially concentrated on 
5 "software development" or building of the formal theory of quantum computing. 
Software development for quantum computing involves attempting to set the 
Hamiltonian of a quantum system to correspond to a problem requiring solution. 
Milestones in these efforts were the discoveries of the Shor and Grover algorithms. 
See, e.g., P. Shor, SI AM J. ofComput., 26:5, 1484-1509 (1997); L. Grover, Proc. 

1 0 28th STOC, 212-219 (1996); and A. Kitaev, LANL preprint quant-ph/951 1026 

(1995). In particular, the Shor algorithm permits a quantum computer to factorize 
natural numbers. Showing that fault-tolerant quantum computation is theoretically 
possible opened the way for attempts at practical realizations of quantum 
computers. See, e.g., E. Knill, R. Laflamme, and W. Zurek, Science, 279, p. 342 

15 (1998). 

One proposed application of a quantum computer is factoring of large 
numbers. In such an application, a quantum computer could render obsolete all 
existing encryption schemes that use the "public key" method. In another 
application, quantum computers (or even a smaller scale device such as a quantum 
2 0 repeater) could enable absolutely safe communication channels where a message, 
in principle, cannot be intercepted without being destroyed in the process. See, 
e.g., H.J. Briegel et al., LANL preprint quant-ph/9803056 (1998) and the 
references therein. 

Quantum computing generally involves initializing the states of N qubits 

2 5 (quantum bits), creating controlled entanglements among the N qubits, allowing 

the quantum states of the qubit quantum system to evolve under the influence of 
the entanglements, and reading the qubits after they have evolved. A qubit 
quantum system is conventionally a system having two degenerate quantum states, 
where the state of the qubit quantum system can have non-zero probability of being 

3 0 found in either degenerate state. Thus, N qubit quantum systems can define an 

initial state that is a combination of 2 N states. The entanglements between qubits 


-2- 


and the interactions between the qubits and external influences control the 
evolution of the distinguishable quantum states and define calculations that the 
evolution of the quantum states perform. This evolution, in effect, can perform 2 N 
simultaneous calculations. Reading the qubits after evolution is complete 
determines the states of the qubit quantum systems and the results of the 
calculations. 

Several physical systems have been proposed for the qubits in a quantum 
computer. One system uses chemicals having degenerate nuclear spin states, see 
U.S. Patent No. 5,917,322, "Method and Apparatus for Quantum Information 
Processing", to N. Gershenfeld and I. Chuang. Nuclear magnetic resonance 
(NMR) techniques can read the spin states. These systems have successfully 
implemented a search algorithm, see, e.g., J. A. Jones, M. Mosca, and R. H. 
Hansen "Implementation of a Quantum Search Algorithm on a Quantum 
Computer," Nature, 393,344-346 (1998) and the references therein, and a number 
ordering algorithm, see, e.g., Lieven M. K. Vandersypen, Matthias Steffen, 
Gregory Breyta, Costantino S. Yannoni, Richard Cleve and Isaac L. Chuang, 
"Experimental Realization of Order-Finding with a Quantum Computer," LANL 
preprint quant-ph/0007017 (2000), Phys. Rev. Lett, Vol. 85, No. 25, 5452-55 
(2000) and the references therein. The number ordering algorithm is related to the 
quantum Fourier transform, an essential element of both Shor's algorithm for 
factoring of a natural number and Grover's Search Algorithm for searching 
unsorted databases, see T. F. Havel, S. S. Somaroo, C.-H. Tseng, and D. G. Cory, 
"Principles and Demonstrations of Quantum Information Processing by NMR 
Spectroscopy, 2000," LANL preprint quant-ph/98 12086 V2 (1999), and the 
references therein. However, efforts to expand such systems to a commercially 
useful number of qubits face difficult challenges. 

Another physical system for implementing a qubit includes a 
superconducting reservoir, a superconducting island, and a dirty Josephson 
junction that can transmit a Cooper pair (of electrons) from the reservoir into the 
island. The island has two degenerate states. One state is electrically neutral, but 
the other state has an extra Cooper pair on the island. A problem with this system 


is that the charge of the island in the state having the extra Cooper pair causes long 
range electric interactions that interfere with the coherence of the state of the qubit. 
The electric interactions can force the island into a state that definitely has or lacks 
an extra Cooper pair. Accordingly, the electric interactions can end the evolution 
of the state before calculations are complete or qubits are read. This phenomenon 
is commonly referred to as collapsing the wavefunction, loss of coherence, or 
decoherence. See Y. Nakamura, Yu. A. Pashkin and J. S. Tsai "Coherent Control 
of Macroscopic Quantum States in a Single-Cooper-Pair Box," Nature V. 398 No. 
6730, P.786 - 788 (1999), and the references therein. 

Another physical system for implementing a qubit includes a radio 
frequency superconducting quantum interference device (RF-SQUID). See J. E. 
Mooij, T. P. Orlando, L. Levitov, Lin Tian, Caspar H. van der Wal, and Seth 
Lloyd, "Josephson Persistent-Current Qubit," Science 285, 1036-39 (August 13, 
1999), and the references therein. The energy levels of the RF-SQUID correspond 
to differing amounts of magnetic flux threading the SQUID ring. Application of a 
static magnetic field normal to the SQUID ring may bring two of these energy 
levels, corresponding to different magnetic fluxes threading the ring, into 
resonance. Typically, external AC magnetic fields are also applied to pump the 
system into excited states so as to maximize the tunneling frequency between qubit 
basis states. A problem with this system is that the basis states used are not 
naturally degenerate and the required biasing field has to be extremely precise. 
This biasing is possible for one qubit, but with several qubits, this bias field fine- 
tuning becomes extremely difficult. Another problem is that the basis states used 
are typically not the ground states of the system, but higher energy states populated 
by external pumping. This requires the addition of an AC field generating device, 
whose frequency will differ for each qubit as the individual qubit parameters vary. 

The race to create the first scalable, practical, and powerful solid state 
quantum computer has existed for over ten years. Ever since the notion of a 
quantum computer first became evident with Feynman in 1982, scientists have 
been creating qubits of various forms. There are currently a number of disclosed 
qubits, where the quantum states are realized in the doubly degenerate ground 


states of the flux in a superconducting loop. Inevitably, these qubits are only 
useful when controlled by magnetic fields, or by some other means which couple 
the qubit to the environment or provide other potential sources of decoherence. In 
order to overcome these sources of decoherence, a large amount of overhead is 
required to control and harvest the quantum power available from the qubit. 
However, the means by which this can be accomplished has as yet eluded 
scientists. Thus, there is a need for a qubit which does not require the coupling 
magnetic fields, but which can be controlled by applying and reading currents and 
voltages. 

There therefore exists a need for integrated solid state structures that can 
form the basic building blocks out of which integrated circuits using quantum 
effects can be built. The desired structures are such that they can be read from, 
written to and operated on in an efficient and scalable manner. 

SUMMARY 

In accordance with the present invention, a qubit is comprised of a multi- 
terminal junction, where two of the terminals of the junction are directly connected 
together, thus forming a superconducting loop. The superconducting loop 
introduces a phase shift so that the phase of the superconducting order parameter *F 
is shifted by an in transition through the structure, where a ranges from -1, 
through zero (no phase shift), to 1. A phase shift can be produced, for example, by 
the inclusion of a phase shifter in the superconducting loop or by external 
application of a magnetic field through the superconducting loop. 

A qubit according to the present invention can be constructed from a multi- 
terminal Josephson junction in which at least two terminals of the junction are 
coupled to a superconducting loop to form a superconducting loop and at least two 
further terminals are open and can be coupled to external current sources. The 
multi-terminal junction can be made of superconducting leads coupled, for 
example, by any combination of constriction junctions (also referred to as micro- 
bridges), tunnel junctions, or semiconducting two dimensional electron gas 


structures within the physical location. In some embodiments of the qubit, the 
terminals of the multi-terminal junction are coupled in a physical location whose 
size is less than the size of the qubit. 

In some embodiments of the invention, properties of both a symmetric 
junction and an asymmetric junction can be utilized. In a symmetric junction, a 
change in the direction of the transport current in the junction equally affects 
current in the terminals of the superconducting loop, thus having no overall affect 
on the current in the loop. In an asymmetric junction, a change in the direction of 
the transport current differentially affects the terminals that form the 
superconducting loop, thus changing the overall current in the loop. 

A symmetric junction allows for the reduction of the potential energy 
barrier between the two nearly degenerate ground states of the quantum system of 
the qubit, thus providing a means of applying a cr x quantum gate operation. An 
asymmetric junction allows for biasing of one of the two ground states of the qubit, 
thus providing a means of applying a tr z quantum gate operation. 

A phase shifter is any structure that shifts the phase of the superconducting 
order parameter *F by an in transition through the structure, where a is a constant 
such that -1 < a < 1 . The phase shift in the superconducting loop causes time- 
reversal symmetry breakdown in the qubit quantum system and thus causes a 
double degeneracy of the ground state without requiring an external magnetic flux 
or other influence. In some embodiments, the terminals in a multi-terminal 
junction can be physically asymmetric. This asymmetry affects the properties of a 
qubit according to the present invention by controlling the phase shift of the order 
parameter *F in transition through a multi-terminal junction. 

A qubit according to the present invention may be constructed out of any 
superconducting material. Embodiments of qubits having any desired number of 
terminals and a phase shifter can also be constructed in accordance with desired 
applications for the qubit. Embodiments of qubit structures include, for example, 
s-wave superconductor / normal metal / d-wave superconductor / normal metal / s- 
wave superconductor, referred to as S-N-D-N-S junctions, superconductor / 
ferromagnet / superconductor, referred to as S-F-S junctions, s-wave 


superconductor / two dimensional electron gas / s-wave superconductor, referred to 
as S-2DEG-S junctions, or multi-crystal d-wave superconductors patterned on an 
insulating substrate. The equilibrium ground state of the qubit quantum system is, 
in the absence of external magnetic fields, twice degenerate, with one of the energy 
levels corresponding to a magnetic flux threading the loop in one sense 
(corresponding to an equilibrium supercurrent flow, for example, in the clockwise 
direction around the superconducting loop), and the other energy level 
corresponding to a magnetic flux threading the loop in the opposite sense 
(corresponding to an equilibrium supercurrent flow, for example, in the 
counterclockwise direction around the superconducting loop). 

Some embodiments of qubits according to the present invention include an 
s-wave (for example, niobium, aluminum, lead, mercury, or tin) superconducting 
structure that includes an asymmetric four-terminal junction with all terminals 
connected by constriction junctions. Two of the terminals can be joined to form a 
superconducting loop and the other two terminals can be coupled to a source of 
transport current. The superconducting loop includes a phase shifter, which may 
consist of a S-N-D-N-S (for example, niobium / gold / Y Ba 2 Cu 3 0 7 . x / gold / 
niobium) junction. If the incoming current is parallel to the a (or b) 
crystallographic direction of the d-wave material, and the outgoing current is 
parallel to the b (or a) crystallographic direction of the d-wave material, this S-N- 
D-N-S junction can give a phase shift of tt . Choosing the incoming and outgoing 
currents to be at any arbitrary angle to each other in the a-b plane in this 
embodiment allows a more general phase shift. 

A magnetic field may also be applied to the superconducting loop. Both 
the transport current and the external magnetic field may be controlled so as to 
initialize the state of the qubit quantum system, allow control of the evolution of 
the qubit quantum system state and read the final state of the qubit quantum system 
after the evolution (and therefore the desired calculation) is complete. Further, 
qubits can be selectively entangled by coupling superconducting loops from 
different qubit structures with a switchable junction, allowing for control of 
entanglements in a qubit array. 


A qubit according to the present invention can include a junction with any 
number of terminals. Some embodiments of the invention include a five terminal 
junction. A superconducting loop is formed between two terminals of the five 
terminal junction. The remaining three terminals, two terminals adjacent to the 
5 looping terminals, and one terminal centrally opposite the looping terminals, form 
a means by which to implement all desired quantum operations, including reading, 
writing, a cr x gate operation, and a a z gate operation. 

Some embodiments, such as the five-terminal qubit, include both 
symmetric and asymmetric properties. Because the critical current in the junction 

1 0 depends on the state of the qubit, a read operation can be performed by applying a 
current asymmetrically across the junction, with a magnitude between the critical 
currents of the two states and determining if a resistance is created. Additionally, a 
ct z gate operation can be performed by applying a pulse of current asymmetrically 
across the junction while a cj x gate operation can be performed by applying a pulse 

15 of current symmetrically across the junction. 

In accordance with some embodiments of the invention, a quantum 
computing method cools a structure containing at least one multi-terminal qubit to 
a temperature that makes the structure superconducting and suppresses 
decoherence processes in the system. The actual temperature will depend on the 

2 0 superconducting materials of the qubit. After the structure is at the appropriate 
temperature, a supercurrent can be established in each superconducting loop, the 
supercurrent being in a particular classical bit state corresponding to the 
information input for the quantum calculation. The quantum systems of each of 
the plurality of qubits is then allowed to evolve in the presence of externally 

2 5 applied magnetic fields and transport currents (whose details constitute the 

"software" or algorithm being followed by the structure). This allows each 
superconducting loop bit state to evolve into quantum states that are admixtures of 
a first state having a first magnetic moment and a second state having a second 
magnetic moment. These quantum states evolve under the action of the system's 

3 0 Hamiltonian in the manner prescribed by quantum mechanics. The evolution 

performs the quantum computation. Determining a measured magnetic moment or 


flux due to the supercurrent in each superconducting loop determines the result of 
the quantum computation. 

In accordance with another aspect of the invention, determining the 
measured magnetic moments of the quantum state on the qubit can also include 
applying an alternating transport current and/or magnetic bias field to each qubit 
and then measuring the magnetic flux produced by the supercurrent flowing in 
each superconducting loop. In some embodiments, a static transport current and/or 
magnetic bias field can be applied to each qubit and the voltage across at least two 
of the terminals in the multi-terminal junction measured to determine the quantum 
state of the qubit. In some embodiments, the quantum states of the qubits can be 
read directly with, for example, a SQUID magnetometer. 

In further aspects of the invention, quantum qubits can be selectively 
entangled with switchable junctions. In embodiments where qubits include a 
superconducting loop, an array of qubits can be entangled by switchably coupling 
the superconducting loops of the array. Additionally, a switchable junction can be 
included to decouple selected ones of the superconducting loops from other multi- 
terminal junctions. 

These and other embodiments according to the present invention are further 
discussed below with respect to the following figures. 

BRIEF DESCRIPTION OF THE FIGURES 

Figures 1A, IB, ID, and IF through 1M show plan views of embodiments 
of multi-terminal junctions in accordance with embodiments of the present 
invention. 

Figures 1C and IE show cross-sectional views of embodiments of multi- 
terminal junctions in accordance with embodiments of aspects of the present 
invention. 

Figures 2A-2G show example embodiments of phase shifter structures in 
accordance with aspects of the present invention. 


Figure 3A shows a plan view of a four-terminal qubit with intrinsic phase 
shifter qubit horizontal architecture in accordance with embodiments of the present 
invention. 

Figure 3B shows a plan view of a multi-terminal qubit with intrinsic phase 
5 shifter qubit horizontal architecture in accordance with embodiments of the present 
invention. 

Figure 3C shows a plan view of a plurality of multi-terminal junctions 
connected in a superconducting loop with intrinsic phase shifter qubit horizontal 
architecture in accordance with embodiments of the present invention. 
1 0 Figure 4A shows a plan view of a plurality of multi-terminal qubits, each 

having a superconducting loop with an intrinsic phase shifter in accordance with 
embodiments of the present invention. 

Figure 4B shows a plan view of a two-terminal circuit with a plurality of 
multi-terminal qubits each having a superconducting loop with intrinsic phase 
15 shifter in accordance with embodiments of the present invention. 

Figure 5 shows a plan view of a six -terminal junction connecting two 
superconducting loops, each with an intrinsic phase shifter, to form a pair of 
qubits in accordance with embodiments of the present invention. 

Figures 6 A and 6B show plan views of a plurality of multi-terminal 
2 0 junctions coupling a plurality of superconducting loops with intrinsic phase shifters 
to form a plurality of qubits in accordance with embodiments of the present 
invention. 

Figure 7 shows a plan view of a voltage measurement circuit in accordance 
with embodiments of the present invention. 

2 5 Figure 8 shows a plan view of a five-terminal quantum qubit according to 

embodiments of the present invention. 

Figure 9 shows an exemplary method of entangling a pair of five terminal 
qubits as shown in Figure 8. 

Figures 10a and 10b illustrate a switchable entanglement method for 

3 0 coupling two five terminal qubits. 
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Figure 1 1 shows an array of N five terminal qubits according to the present 
invention. 

Figure 12 shows another embodiment of a qubit array according to the 
present invention. 

5 Figure 13 shows a plan view of a quantum qubit according to the present 

invention. 

Use of the same reference symbols in different figures indicates elements 
having similar or identical functions. 

10 DETAILED DESCRIPTION 

In accordance with embodiments of the invention, a quantum computing 
operation can be performed on an array of quantum qubits where at least one of the 
qubits includes a qubit according to the present invention. A qubit according to the 
present invention includes a multi-terminal junction where two terminals of the 
1 5 multi-terminal junction are joined to form a superconducting loop. The 

superconducting loop introduces a phase shift to the superconducting order 
parameter. In some embodiments, the superconducting loop includes an intrinsic 
phase shifter. 

Intrinsic phase shifters in superconducting phase quantum bits (qubits) are 

2 0 disclosed in M.H.S. Amin, T. Duty, A. Omelyanchouk, G. Rose and A. Zagoskin, 

U.S. Provisional Application Serial No. 60/257624, "Intrinsic Phase Shifter as an 
Element of a Superconducting Phase Quantum Bit", filed December 22, 2000, 
herein incorporated by reference in its entirety. A phase shifting structure with 0 
and n -phase shifts in a two-terminal DC SQUID is described in R.R. Schulz, B. 
25 Chesca, B. Goetz, C. W. Schneider, A. Schmehl, H. Bielefeldt, H. Hilgenkamp, J. 
Mannhart and C. C. Tsuei, "Design and Realization of an all d-Wave dc n- 
Superconducting Quantum Interference Device", Appl. Phys. Lett. 76, 7 p.912-14 
(2000), herein incorporated by reference in its entirety. 

Figure 13 shows a block diagram of a qubit 100 according to the present 

3 0 invention. Qubit 100 includes multi-terminal junction 120 with terminals 110-1 

through 1 10-N where N is an integer. At least two of terminals 110-1 through 110- 
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N are coupled to portions 124 and 125 of a superconducting loop 122. In Figure 
13, terminals 110-1 and 1 10-N are coupled to superconducting loop 122, but in 
general any two of terminals 1 10-1 through 1 10-N can be coupled to 
superconducting loop 122. Superconducting loop 122 further includes a phase 
shifting structure 123 coupled into superconducting loop 122. In Figure 13, 
terminals 110-i and 110-j refer to arbitrary ones of terminals 110-1 through 110-N. 
Although terminals 110-1 and 1 10-N are shown coupled to portions 124 and 125, 
respectively, of superconducting loop 122, in general any pair of terminals 1 10-1 
through 1 10-N can be coupled to terminals 124 and 125 of superconducting loop 
122. 

Embodiments of qubits 100 according to the present invention can include 
phase shifter 123, which introduces an arbitrary phase shift inclusively between -n 
and 71. Phase shifters, such as phase shifter 123, that introduce arbitrary phase 
shifts can be more practical for the construction of qubits since the phase shifters 
with arbitrary phase shifts are more easily constructed. Further tuning of the phase 
shift accumulated through superconducting loop 122 in qubit 100 can be 
accomplished by the application of a magnetic field or adjustment of the transport 
current I T in superconducting loop 122. 

Additionally, embodiments of the present invention include at least one 
terminal junction 120. Terminal junction 120 joins at least two terminals, 
terminals 110-1 through 110-N. In some embodiments, the physical size of 
junction 120 is much less than the size of superconduction loop 122. 

Four-terminal SQUID devices are discussed in A. N. Omelyanchouk and 
Malek Zareyan, "Ballistic Four-Terminal Josephson Junction: Bistable States and 
Magnetic Flux Transfer", Los Alamos preprint cond-mat/9905139, and BJ. 
Vleeming, "The Four-Terminal SQUID", Ph.D. Dissertation, Leiden University, 
The Netherlands, 1998, both of which are herein incorporated by reference in their 
entirety. Four terminal SQUID devices are further discussed in R. de Bruyn 
Ouboter and A. N. Omelyanchouk, "Macroscopic Quantum Interference Effects in 
Superconducting Multiterminal Structures", Superlattices and Microstructures, 
Vol. 25 No 5/6 (1999), herein incorporated by reference in its entirety. 


A quantum computation relies on qubit 100 including a qubit quantum 
system formed by supercurrents on superconducting loop 122 having degenerate 
ground states, designated j0> and |1>, of the supercurrent. An array of multi- 
terminal superconducting loops 122 having phase shifters 123 can be fabricated in 
useful numbers in a solid state structure. The ground state of the qubit quantum 
system includes two states that correspond to supercurrent flows that circulate 
clockwise and counterclockwise, respectively, in the plane of superconducting loop 
122. The qubit quantum system of qubit 100 can be initialized by the introduction 
of supercurrents from an external source some of terminals 110-1 through 1 10-N 
not coupled to superconducting loop 122. The ground-state of the qubit quantum 
system in each superconducting loop 122 containing phase shifter 123 is doubly 
degenerate in the absence of externally applied magnetic fields and/or transport 
currents I T (each circulation direction has the same energy) and provides the basis 
for a qubit 100 for quantum computing in accordance with embodiments of the 
present invention. 

The two degenerate states, corresponding to classical bit states, represented 
symbolically as |0> and |1>, are then the two basis states of the qubit quantum 
system of qubit 100. The magnitude of the flux threading superconducting loop 
122 can be much less than half a flux quantum <D 0 , both because of intrinsic phase 
shifter 123 and the presence of the terminals 110-1 through 1 10-N, which also 
introduces a phase shift. At least two external terminals, terminals 110-2 and 1 10-j 
in Figure 12, for example, of qubit 100 can be coupled to sources of transport 
current Ii and Ij, respectively, in Figure 12. Terminals 1 10-i and 110-j along with 
the current source creates a transport current loop 127. Additionally, an external 
magnetic flux, indicated by field B , can be applied through superconducting loop 
122 in order to control the physical parameters of the qubit quantum system of 
qubit 100. By changing the transport current I T and/or applying an external 
magnetic field B , the magnitude of the flux threading superconducting loop 122, 
the potential barrier between the two basis states |0> and |1> of the qubit quantum 
system of qubit 100, and the tunneling matrix element A T (I) between the basis 
states of the qubit quantum system can be adjusted. 
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The choice of the physical sizes of the constriction junctions of junction 
120, tunnel junctions and/or semiconducting two-dimensional electron gas 
structures that couple the terminals, also affects the functioning of qubit 100. To 
achieve a small total flux in superconducting loop 122 (which is desirable for 
5 decreasing the decoherence rate) and maximum influence of the transport current 
It on the properties of superconducting loop 122, in some embodiments the links in 
the transport loop (e.g., the current loop providing current to junction 120) are 
much wider than the ones in superconducting loop 122. A small residual flux 
exists because of spontaneous supercurrents, even in the absence of external fields. 

10 In those embodiments, the height of the potential energy barrier between the two 
degenerate quantum states of the qubit quantum system of qubit 100 will be 
affected most pronouncedly by the transport currents I T . 

Multi-terminal junction 120 includes two important regimes: symmetric 
and asymmetric. Further, multi-terminal junction 120 can display symmetric, 

1 5 asymmetric, or a combination of symmetric and asymmetric properties. In a 
symmetric junction, a change in the direction of the transport current in the 
junction equally affects current in the terminals of the superconducting loop, thus 
having no overall affect on the current in the loop. Whereas, in an asymmetric 
junction, a change in the direction of the transport current differentially affects the 

2 0 terminals that form the superconducting loop, thus changing the overall current in 
the loop. 

A symmetric junction can be used to reduce the potential energy barrier 
between the two nearly degenerate ground states of qubit 100, thus providing a 
means of applying a ct x quantum gate operation. If a change in the direction of the 

2 5 transport current in terminals 1 1 0-2 through 1 1 0-j has an equal effect on the 

current in terminals 110-1 and 1 10-N, then junction 120 is symmetric with respect 
to terminals 110-2 and 1 10-j. Qubit 100, under these circumstances, can then be 
referred to as a "symmetric qubit". 

An asymmetric junction can be used to bias one of the two ground states of 

3 0 the qubit, thus providing a means of applying a o z quantum gate operation. If a 

change in the direction of the transport current in terminals 110-2 through 1 10-j 
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causes a differential change in the current in the loop terminals 1 10-1 and 1 10-N, 
then the junction is said to be asymmetric. Qubit 1 00, then, can be referred to as 
an "asymmetric qubit". 

Therefore, a four-terminal junction 120 can be either a symmetric or 
asymmetric junction, and four-terminal qubit 100 can then be either a symmetric or 
asymmetric qubit. If qubit 100 includes a junction 120 with more than four 
terminals, for example a five-terminal qubit, then both symmetric and asymmetric 
properties can be realized. 

An asymmetric qubit can be written to (i.e., the quantum states initialized) 
by applying a transport current I T when the magnitude of the transport current I T is 
larger than a threshold value which is determined by the specific implementation 
of qubit 100. The direction of transport current I T is chosen depending on which 
basis state (i.e., |0> or jl>) is being written into the qubit quantum system. The 
application of transport current I T , then, has the effect of biasing the qubit quantum 
system states into one of the degenerate basis states. In the biased state the qubit 
quantum system will decay to the most energetically favorable state (either |0> or 
|1> as required). In such systems, the time to decay typically is shorter than about 
1 millisecond, depending on the particular embodiment of qubit 100. Depending 
on the particular embodiment of qubit 100, a magnetic field B can also be applied 
in one of two directions, which can alter the time to decay. The magnetic field can 
be applied opposing the transport current induced bias, thereby decreasing the time 
for decay, or supporting the transport current induced bias, thereby increasing the 
time to decay. 

A symmetric qubit can be written to (i.e., biased) by applying a static 
magnetic field B . As described above, this will cause the qubit quantum system 
of qubit 100 to decay into the energetically favorable state on a time-scale 
dependent upon the embodiment of qubit 100 and the magnitude of an externally 
applied magnetic field B . 

Single qubit operations on asymmetric qubits can be performed by 
modulating the transport current and/or the external magnetic field strength. 
Setting the transport current I T to zero sets the effective Hamiltonian describing the 
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quantum system of qubit 100 proportional to d x , which is referred to as a Pauli 
matrix. In the basis where the qubit basis states |0> and |1> are chosen so that the 
state |0> corresponds to the vector (1,0) and the state |1> corresponds to the vector 

(0,1), cf x = I I . This basis can be called the Z-diagonal basis. In this basis the 

5 Pauli matrix a x rotates one of the basis states into the other basis state (i.e., 
<? x |0>=|1> and a x |1>=|0>). 

Increasing the transport current I T past a threshold current, which can be an 
implementation dependent critical value, sets the Hamiltonian proportional to cf z , 
which is another Pauli matrix. The matrix cf z is defined in the Z-diagonal basis to 

. ~ T 1 °1 

10 be a 2 = I I . The Pauli matrix cf z biases the states, or in an alternative 

interpretation adds a phase to the second state (i.e. cf z |0>=|0> and cf z |1>=-|1>). 
An arbitrary single qubit operation can be performed by performing combinations 
of the functions described by cf x and cf z . 

To keep the quantum system of qubit 100 in some specific state an 

15 alternating transport current I T (?) can be applied. The Hamiltonian representing 
the quantum system of qubit 100, then, is proportional toI T (t)a z . In some 
embodiments, for example, I T (t) can be a square wave. This method can be used 
in conjunction with a clock whose frequency is an integer multiple of the 
frequency of I T (t) so that, at every clock pulse, the quantum system of qubit 100 

2 0 is in the same state in which it began. In addition, during the evolution of the qubit 
states external magnetic fields B maybe applied in accordance with a specific 
usage of the qubit. 

Single qubit operations on symmetric embodiments of qubit 100 can be 
performed by modulating the transport current I T (t) and/or the external magnetic 

25 field B . In the symmetric qubit embodiments of qubit 100, changing the transport 
current I T (?) does not affect the bias (i.e., the initial quantum states of the qubit 
quantum system) but does affect the energy barrier between the qubit quantum 


system basis states. The effective Hamiltonian describing the quantum system of 
qubit 100, then, includes a term proportional to A T (l)a x where the tunneling 
matrix element A r (/) can be varied over a large range dependent on the transport 
current I T (*) . Applying a magnetic field B normal to the plane of the 
superconducting loop 122 provides another term of the Hamiltonian that is 
proportional to 6 2 . To keep the quantum system in some specific state, an 
alternating magnetic field normal to the superconducting loop 122, B(t) , can be 
applied, adding a term proportional to B(t)6 z to the Hamiltonian. In some 
embodiments, for example, B(t) can be modulated in a square wave. This method 
can also be used in conjunction with a clock whose frequency is an integer 
multiple of the frequency of B(t) so that at every clock pulse the qubit quantum 
system of qubit 100 is in the same state in which it began. 

In some embodiments, the state of the qubit quantum system of an 
asymmetric embodiment of qubit 100 begins with the application of an alternating 
transport current I T (t) as described above. The magnetic flux threading the 
superconducting loop 122 can then be measured using a suitable magnetic field 
measuring device, for example a SQUID magnetometer or a magnetic force 
microscope. The alternating transport current I T (t) causes the decay rate out of 
the desired final state to be minimized. Once the state of the qubit quantum system 
of qubit 100 is measured, the transport current I T (t) can be set to bias the states 
such that the measured state has a lower energy level, making the measured state's 
information stable so that it may be accessed at any later time. 

In some embodiments, the state of the qubit quantum system of the 
asymmetric qubit can be measured by applying a transport current I T {t) in a fixed 
direction and then monitoring the voltage drop across pairs of external leads, for 
example 110-2, 1 10-i and 1 10-j of Figure 13. If the measured voltage remains 
constant then the state of the qubit quantum system corresponds to the favored 
energy state indicated by the fixed direction of the transport current I T (t) . If the 
measured voltage changes, then the state of the qubit quantum system corresponds 


to an excited state indicated by the biasing of the fixed direction of the transport 
current I T (t) . The changing voltage indicates a decay of the qubit quantum 
system state from the excited state. This indicates that the qubit quantum system 
was in the excited state relative to the bias indicated by the fixed direction of the 
transport current at the end of the calculation. 

In some embodiments, the state of the quantum system of a symmetric 
qubit can be read similarly to the state of the quantum system of an asymmetric 
qubit, except that biasing is performed by application of magnetic fields normal to 
superconducting loop 122 and not via application of transport currents. 

Figure 1 A shows a plan view of an embodiment of a four-terminal 
constriction junction 120 according to the present invention. Four-terminal 
junction 120 includes terminals 110-1, 110-2, 1 10-3, and 110-4 coupled at a 
constriction junction 140. Superconducting currents Ii, I 2 , 13 and I4 can exist in 
terminals 110-1, 110-2, 110-3 and 110-4, respectively. The terminal lengths Li 
through L 8 , which describe the linear dimensions of terminals 110-1, 110-2, 1 10-3, 
and 1 10-4 as indicated in Figure 1 A, can all be different and are typically chosen to 
be less than about 10 microns. The terminal widths Wi through W 4 , which 
describe the widths of terminals 110-1 through 1 10-4, respectively, can also all be 
different and are typically chosen to be less than the coherence length of the 
superconducting material of terminals 110-1 through 1 10-4. For example, if 
terminals 1 10-1 through 1 10-4 are of aluminum, the coherence length is 1.6 
microns. Four-terminal constriction junction 120 can be fabricated of any 
superconducting material. 

In an exemplary embodiment, four-terminal constriction junction 120 can 
be fabricated of aluminum. Widths Wi and W 2A can each be approximately 0.5 
microns; widths W 3 and W 4 can each be approximately 0.05 microns; lengths Li, 
L 4 , L 5 , L 6 , L 7 , and L 8 can each be approximately 1 micron; and lengths L 2 and L 3 
can each be approximately 0.55 microns. 

In another exemplary embodiment, four-terminal constriction junction 120 
can be fabricated of aluminum, with width Wi approximately 0.5 microns, W 2 
approximately 0.3 microns, W 3 approximately 0.08 microns, and W 4 


approximately 0.05 microns, lengths Li, L 4 , L 5 , L 6 , and L 7 each approximately 1 
micron, L 2 approximately 0.75 microns, L 3 approximately 0.68 microns, and L 8 
approximately 0.9 microns. 

Figure IB shows a plan view of another embodiment of a four-terminal 
5 junction 120 with a constriction junction 142 coupling terminals 1 10-3 and 1 10-4 
and two tunnel junctions 140 and 141. Four-terminal junction 120 of Figure IB 
further includes terminals 110-1 and 1 10-2. Tunnel junction 140 couples terminal 
1 10-4 with terminal 1 10-1. Tunnel junction 141 couples terminal 110-3 with 
terminal 110-2. Superconducting currents Ii, I 2 , 1 3 and I4 can exist in terminals 

10 110-1, 110-2, 110-3 and 110-4, respectively. The terminal lengths Li - L 9 , which 
describe the linear dimensions of terminals 110-1 through 110-4 and the separation 
between junctions 140 and 141, can all be different and are typically chosen to be 
less than about 10 microns. The terminal widths Wj - W 4 , which as before 
describe the widths of terminals 110-1 through 110-4, respectively, can also be 

1 5 different and are typically chosen to be less than the coherence length of the 

superconductor used. For example, in the case of aluminum, the coherence length 
is about 1.6 microns. Four-terminal junction 120 can be fabricated of any 
superconductor material. Tunnel junction 140 is typically fabricated using an 
insulating layer between terminals 110-4 and 110-1. Tunnel junction 141 is 

2 0 typically fabricated using an insulating layer between terminals 1 1 0-2 and 1 1 0-3 . 

Figure 1C shows a cross-sectional diagram of the embodiment of four- 
terminal junction 120 shown in Figure IB. Terminals 110-4 and 110-1 are 
separated by insulator 150 at junction 140 and terminals 110-3 and 110-2 are 
separated by insulator 150 at junction 141. 
25 In an exemplary embodiment, four-terminal junction 120 as shown in 

Figure IB can be fabricated of aluminum, with widths Wi and W 2 being 
approximately 0.05 microns, widths W 3 and W 4 being approximately 0.5 microns, 
lengths Li, L 2 , L 3 , L 4 , L 5 , and L 7 being approximately 1 micron, and L 9 being 
approximately 0.1 microns. 

3 0 Tunnel junction 140 can be fabricated where insulation layer 1 50 (Figure 

1C) is an aluminum oxide layer of thickness approximately 0.05 microns. 
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Insulating layer 150 can be deposited onto terminal 110-4 Terminal 1 10-1 is then 
deposited on insulating layer 150. The tunnel junction 141 can also be fabricated 
where insulating layer 150 (Figure 1C) is an aluminum oxide layer of thickness 
approximately 0.05 microns which is deposited onto terminal 110-3, upon which is 
5 deposited terminal 110-2. Although in Figure 1C, terminals 1 10-1 and 1 10-2 are 
shown above terminals 110-4 and 110-3, respectively, in general terminals 110-1 
and 110-2 can be on either side of terminals 1 10-4 and 110-3, respectively. 

In another exemplary embodiment, four-terminal junction 120 as shown in 
Figures IB and 1C can be fabricated of aluminum with the following dimensions: 

1 0 width Wi being approximately 0.05 microns; width W 2 being approximately 0.08 
microns; width W 3 being approximately 0.3 microns; width W 4 being 
approximately 0.5 microns; lengths L 2 , L 3 , L 4 , L 5 , andL 7 being approximately 1 
micron; length Li being approximately 0.8 microns; and length L 9 being 
approximately 0.1 microns. Again, junctions 140 and 141 can be fabricated with 

1 5 insulation layer 150 being an approximately 0.05 microns layer of aluminum oxide 
deposited onto terminals 1 10-4 and 110-3, respectively. Terminals 1 10-1 and 1 10- 
2 can then be deposited at junctions 140 and 141 on insulating layer 150. 

Figure ID shows a plan view of another embodiment of four-terminal 
junction 120 with one constriction junction 142 coupling terminals 110-3 and 110- 

2 0 4 and tunnel junctions 140 and 141 , which are fabricated so as to run parallel to 

terminals 1 10-4 and 110-3 over distances Di and D 2 , respectively. In some 
embodiments, Di and D 2 are each greater than approximately 0.2 microns. Tunnel 
junction 140 couples terminal 110-4 with terminal 110-1. Tunnel junction 141 
couples terminal 1 10-2 with terminal 110-3. Superconducting currents Ii, I 2 , 1 3 and 
25 I4 are the currents in terminals 110-1,11 0-2, 1 1 0-3 and 1 1 0-4, respectively. The 
terminal lengths L 2 , L4, L5 and L7 can all be different and are typically chosen to be 
less than about 10 microns. Note that in some embodiments the length L 9 can be 
negative, which may result in overlapping tunnel junctions 140 and 141. 
Terminals 110-1 and 110-2 can extend at arbitrary angles Ai 20 -i and A J2 o- 2 

3 0 respectively. The terminal widths Wi - W 4 can all be different and are typically 

chosen to be less than the coherence length of the superconducting material used. 
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For example, in the case of aluminum terminals, the coherence length is about 1 .6 
microns. Four-terminal junction 120 can, however, be fabricated of any 
superconducting material. As shown in Figure IE, tunnel junction 140 can be 
fabricated using an insulating layer 150 between terminals 1 10-4 and 1 10-1. 
Tunnel junction 141 can be fabricated using an insulating layer 151 between 
terminals 110-3 and 110-2. 

In an exemplary embodiment, four-terminal junction 120 as shown in 
Figures ID and IE can be fabricated of aluminum with the following dimensions: 
widths Wi and W 2 being approximately 0.05 microns; widths W 3 and W 4 being 
approximately 0.5 microns; lengths L 2 and L 4 being approximately 1 micron; 
lengths Di and D 2 being approximately 1 micron; lengths L 5 and L 7 being 
approximately 1.5 microns; and L 9 being approximately 0.1 microns. Tunnel 
junction 140 can be fabricated with insulating layer 150 being an aluminum oxide 
layer of thickness approximately 0.05 microns deposited onto terminal 1 10-4, upon 
which is deposited terminal 1 10-1. Tunnel junction 141 can be fabricated with 
insulating layer 151 being an aluminum oxide layer of thickness approximately 
0.05 microns which is deposited onto terminal 110-3, upon which is deposited 
terminal 110-2. 

In another exemplary embodiment of junction 120 as shown in Figures ID 
and IE, four-terminal junction can be fabricated of aluminum with the following 
dimensions: width Wi being approximately 0.05 microns; width W 2 being 
approximately 0.08 microns; width W 3 being approximately 0.3 microns; width W 4 
being approximately 0.5 microns; lengths L 2 and L 4 being approximately 1 micron; 
lengths Dl and D2 being approximately 1 micron; lengths L 5 and L 7 being 
approximately 1.5 microns; and length L 9 being approximately 0.1 microns. 
Again, insulating layer 150 forming part of tunnel junction 140 can be fabricated 
of aluminum oxide with thickness approximately 0.05 microns and insulating layer 
151 forming part of tunnel junction 141 can be fabricated of aluminum oxide of 
thickness approximately 0.05 microns. 

Figure IF shows a plan view of an embodiment of a four-terminal junction 
120 with constriction junction 141 coupling terminals 110-2, 1 10-3 and 110-4 and 


tunnel junction 140 coupling junctions 1 10-4 with 110-1. Superconducting 
currents Ii, I 2 , 1 3 and I4 can exist in terminals 110-1,1 10-2, 1 1 0-3 and 1 1 0-4, 
respectively. The terminal lengths Li, L 2 , L 3 , L 4 , L 5 and L 7 , which indicate the 
dimensions of terminals 110-1 through 1 10-4, can all be different and are typically 
chosen to be less than about 10 microns. The widths of terminals 110-1 through 
110-4, Wi - W 4 respectively, can also be different and are typically chosen to be 
less than the coherence length of the superconductor used. Junction 120 can be 
fabricated of any superconducting material. Tunnel junction 140 can typically be 
fabricated using an insulating layer 150 between terminals 110-4 and 1 10-1, as was 
shown for junction 140 in Figure 1C. 

In an exemplary embodiment, four-terminal junction 120 as shown in 
Figure IF can be fabricated of aluminum with the following dimensions: widths 
Wi, and W 4 being approximately 0.5 microns; widths W 2 and W 3 being 
approximately 0.05 microns; lengths L h L 2 , L 3 , L 4 , L 5 , and L 7 being approximately 
1 micron; and length L9 being approximately 0.1 microns. Tunnel junction 140 
can be fabricated with insulating layer 150 being an aluminum oxide layer of 
thickness approximately 0.05 microns. 

In another exemplary embodiment of junction 120 as shown in Figure IF, 
four-terminal junction can be fabricated of aluminum with the following 
dimensions: width Wi being approximately 0.05 microns; width W 4 being 
approximately 0.08 microns; width W 2 being approximately 0.3 microns; width W 3 
being approximately 0.5 microns; lengths Li, L 2 , L 3 , L 4 , L 5 , and L 7 being 
approximately 1 micron; and length L 9 being approximately 0.1 microns. Tunnel 
junction 140 can, as before, include insulating layer 150 fabricated using an 
aluminum oxide layer of thickness approximately 0.05 microns. 

Figure 1G shows a plan view of another embodiment of a four-terminal 
junction 120 having constriction junction 141 coupling terminals 110-2, 110-3 and 
1 10-4 and tunnel junction 140, which is fabricated so as to run parallel to terminal 
110-4 at least over distance D i3 coupling terminals 110-4 with 110-1. In some 
embodiments, distance D 2 is greater than approximately 0.2 microns. As before, 
terminals 1 10-1, 110-2, 110-3 and 110-4 can carry superconducting currents Ii, I 2 , 


I 3 and I4, respectively. Terminal lengths L 2 , L 3 , L 4 , L 5 and L 7 can all be different 
and are typically chosen to be less than about 10 microns. The terminal widths Wi 
- W 4 of terminals 110-1 through 110-4, respectively, can all be different and are 
typically chosen to be less than the coherence length of the superconducting 
material used. Four-terminal junction 120 can be fabricated of any 
superconductor. As discussed before, tunnel junction 140 is typically fabricated 
with an insulating layer 150 between terminals 110-4 and 1 10-1, as shown in 
Figure IE. 

In an exemplary embodiment of junction 120 as shown in Figure 1G, four- 
terminal junction can be fabricated of aluminum with the following dimensions: 
widths Wi and W 4 each being approximately 0.05 microns; widths W 2 and W 3 
each being approximately 0.5 microns; lengths L 2 , L 3 , L 4 , and L 5 each being 
approximately 1 micron; length Di being approximately 0.5 microns; length L9 
being approximately 0.1 microns; and length L 7 being approximately 1.5 microns. 
Tunnel junction 140 can include insulating layer 150 as shown in Figure IE which 
can be fabricated using an aluminum oxide layer of thickness approximately 0.05. 

In another exemplary embodiment of junction 120 as shown in Figure 1G, 
four-terminal junction can be fabricated of aluminum with the following 
dimensions: width Wi being approximately 0.05 microns; width W 4 being 
approximately 0.08 microns; width W 2 being approximately 0.3 microns; width W 3 
being approximately 0.5 microns; lengths L 2 , L 3 , L 4 , and L 5 being approximately 1 
micron, length Di being approximately 0.5 microns, length L 9 being approximately 
0.1 microns; and length L 7 being approximately 1.5 microns. Again, tunnel 
junction 140 can include insulating layer 150 which can be fabricated using an 
aluminum oxide layer of thickness approximately 0.05 microns. 

Figure 1H shows a plan view of an embodiment of a four-terminal junction 
120 where the four terminals 110-1, 110-2, 1 10-3, and 1 10-4 are all coupled by 
junction 140, which in this embodiment is a two dimensional electron gas 
structure. Again, terminals 110-1, 110-2, 110-3, and 110-4 can carry 
superconducting currents Ii, I 2 , 1 3 and I 4 , respectively. The terminal lengths Li - 
Lg can all be different and are typically chosen to be less than about 10 microns. 


The terminal widths Wi - W 4 also can all be different and are typically chosen to 
be less than the coherence length of the superconductor used. Superconducting 
terminals 110-1, 110-2, 110-3 and 1 10-4 can be fabricated of any superconducting 
material which can be coupled to a two dimensional electron gas structure of 
5 junction 140. The two dimensional electron gas structure of junction 140 can be 
fabricated of any structure which allows electrons to be confined to a single two 
dimensional plane and allows coupling of these electrons to superconducting 
terminals 110-1, 110-2, 110-3 and 110-4. 

An exemplary embodiment of four-terminal junction 120 as shown in 

1 0 Figure 1H can be fabricated of niobium with the following dimensions: widths Wi 
and W 2 being approximately 0.2 microns; widths W 3 and W 4 being approximately 
0.05 microns; and lengths Li, L 2 , L 3 , L 4 , L 5 , L 6 , L 7 , and L 8 each being 
approximately 0.5 microns. The two dimensional electron gas structure of junction 
140 can be constructed of a semiconducting InAs heterostructure. The details of 

1 5 the fabrication and behavior of two dimensional electron gas junctions are well 

known, see e.g., A. Jacobs, R. Kiimmel, and H. Plehn, "Proximity Effect, Andreev 
Reflections, and Charge Transport in Mesoscopic Superconducting- 
Semiconducting Heterostructures," Los Alamos Preprint cond.-mat/98 10343 v2 
(1998), republished in Superlattices and Microstructures, Vol. 25, Nr. 5/6, 669- 

2 0 68 1 (1 999), which is herein incorporated by reference in its entirety. 

Another exemplary embodiment of four-terminal junction 120 as shown in 
Figure 1H can be fabricated of niobium with the following dimensions: width Wi 
being approximately 0.1 microns; width W 2 being approximately 0.08 microns; 
width W 3 being approximately 0.03 microns; width W 4 being approximately 0.01 

2 5 microns; and lengths L ]3 L 2 , L 3 , L 4 , L 5 , L 6 , L 7 , and L 8 each being approximately 

0.5 microns. The two dimensional electron gas structure of junction 140 can be 
constructed of semiconducting InAs. 

Figure II shows a plan view of another embodiment of a four-terminal 
junction 120, with constriction junction 141 coupling terminals 110-1 and 110-2 

3 0 and junction 140, which can be a two dimensional electron gas structure, coupling 

terminals 110-4 and 110-3 with constriction junction 141. Terminals 110-1, 110-2, 
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1 10-3, and 110-4 can carry superconducting currents Ii, I 2 , h and Lj, respectively. 
As before, all of the terminal lengths (Li, L 2 , L 4 , L 5 , L 6 , L 7 , L 8 and L 9 ) indicating 
the dimensions of terminals 110-1 through 110-4 can all be different and are 
typically chosen to be less than about 10 microns. The widths of terminals 110-1 
5 through 1 10-4, Wi - W 4 , can all be different and are typically chosen to be less 
than the coherence length of the superconducting material utilized in their 
fabrication. Terminals 110-1, 110-2, 110-3 and 110-4 can be fabricated of any 
superconducting material which can be coupled to a two dimensional electron gas 
junction 140. 

10 In an exemplary embodiment, four-terminal junction 120 as shown in 

Figure II can be fabricated of niobium with the following dimensions: widths Wi 
and W 2 being approximately 0.5 microns; widths W 3 and W 4 being approximately 
0.05 microns; lengths Li, L 2 , L 4 , L 5 , L 6 , L 7 , L 8 , each being approximately 0.5 
microns; and length Lg being approximately 0.1 microns. The two dimensional 

1 5 electron gas structure of junction 140 can be formed of InAs. 

In another exemplary embodiment, four-terminal junction 120 can be 
fabricated of niobium with the following dimensions: width Wi being 
approximately 0.05 microns; width W 4 being approximately 0.08 microns; width 
W 2 being approximately 0.1 microns; width W 3 being approximately 0.15 microns; 

2 0 lengths Li, L 2 , L 4j L 5 , L 6 , L 7 , and L 8 being approximately 0.5 microns; and length L 9 
being approximately 0.1 microns. The two dimensional electron gas structure of 
junction 140 can be formed of InAs. 

Figure 1 J shows a plan view of an embodiment of a four-terminal junction 
120 having tunnel junction 141 coupling terminals 1 10-1 and 110-2 and two 

2 5 dimensional electron gas structure junction 1 40 coupling terminals 1 1 0-3 and 110- 

4 with junction 141. Superconducting currents I], I 2 , 13 and I4 can exist in terminals 
110-1, 110-2, 110-3 and 1 10-4, respectively. The terminal lengths indicating the 
dimensions of terminals 110-1 through 110-4, Li, L 2 , L 4 , L 5 , L 6 , L 7 , L 8 , and Di in 
Figure 1J, can all be different and are typically chosen to be less than 10 microns. 

3 0 The terminal widths Wi - W 4 , indicating the widths of terminals 110-1 through 

1 10-4, can all be different and are typically chosen to be less than the coherence 
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length of the superconducting material of terminals 110-1 through 1 10-4. 
Terminals 110-1 through 1 10-4 can be fabricated of any superconducting material 
which can be coupled to the two dimensional electron gas structure of junction 
140. Tunnel junction 141 is typically fabricated by introducing a layer of 
insulating material, such as aluminum oxide, between terminals 110-1 and 1 10-2, 
as was previously discussed with respect to Figures 1C and IE. 

An exemplary embodiment of four-terminal junction 120 can be fabricated 
of niobium with the following dimensions: widths Wi and W 2 being approximately 
0.5 microns; widths W 3 and W 4 being approximately 0.05 microns; lengths Li, L 2 , 
L 4 , L 5 , L 6 , L 7 , and L 8 each being approximately 0.5 microns; and length Di being 
approximately 0.1 microns. As previously discussed, the two dimensional electron 
gas structure of junction 140 can be formed of InAs. Additionally, tunnel junction 
141 can be fabricated using a niobium oxide layer of thickness approximately 0.05 
microns which is deposited onto terminal 110-2, upon which is deposited terminal 
110-1. 

Another exemplary embodiment of four-terminal junction 120 can be 
fabricated of niobium with the following dimensions: width Wi being 
approximately 0.05 microns; width W 4 being approximately 0.08 microns; width 
W 2 being approximately 0.1 microns; width W3 being approximately 0.15 microns; 
lengths Li, L 2 , L 4 , L5, L6, L 7 , and L 8 each being approximately 0.5 microns; and 
length Di being approximately 0.1 microns. The two dimensional electron gas 
structure of junction 140 can be formed of InAs and tunnel junction 141 can be 
formed with junctions 110-1 and 110-2 separated by a niobium oxide layer of 
approximately 0.05 microns thickness. 

Figure IK shows a plan view of an embodiment of a multi-terminal 
junction 120 with constriction junction 140 coupling terminals 110-1 through 110- 
N, where N is an arbitrary integer. Terminals 110-1 through 1 10-N can carry 
superconducting currents Ii through I N , respectively. The terminal lengths Li 
through L 2N , which describe the dimensions of terminals 110-1 through 110-N, can 
all be different and are typically chosen to be less than about 10 microns. The 
widths of terminals 110-1 through 110-N, Wi - W N , can also all be different and 


are typically chosen to be less than the coherence length of the superconducting 
material of multi-terminal junction 120. Multi-terminal constriction junction 120 
can be fabricated of any superconducting material. In accordance with aspects of 
the present invention, multi-terminal junction 120 can be an element in 
superconducting circuits. 

Figure 1L shows a plan view of an embodiment of a multi-terminal 
junction 120 having a constriction junction 140-2 coupling terminals 1 10-N and 
110-(N-1), a tunnel junction 140-1 coupling terminals 110-N with 110-1, and 
tunnel junctions 140-3 through 140-(N-1) coupling junctions 110-2 through 110- 
(N-2), respectively, to terminal 1 10-(N-1), where N is an arbitrary number. 
Terminals 110-1 through 1 10-N can carry superconducting currents Ii through I N; 
respectively. The lengths L] through L 2N , which characterize the dimensions of 
terminals 110-1 through 1 10-N, can all be different and are typically chosen to be 
less than about 10 microns. The widths of terminals 110-1 through 110-N, Wi - 
W N , can all be different and are typically chosen to be less than the coherence 
length of the superconductor used. The separation between junctions, Di through 
D N . 2 , can also all be different and are typically chosen to be less than about 10 
microns. Multi-terminal junction 120 can be fabricated of any superconducting 
material. In accordance with certain aspects of the present invention, multi- 
terminal junction 120 can be utilized as a circuit element in superconducting 
circuits. 

Figure 1M shows a plan view of another embodiment of multi-terminal 
junction 120. Multi-terminal junction 120 of Figure 1M includes a two 
dimensional electron gas junction 140 coupling terminals 110-1 through 110-N, 
where N is an integer. As before, terminals 110-1 through 1 10-N can carry 
superconducting currents Ii through I N , respectively. The terminal lengths Li 
through L ( 2N), which indicate the dimensions of terminals 110-1 through 1 10-N, 
can all be different and are typically chosen to be less than about 1 0 microns. The 
widths of terminals 1 10-1 through 110-N, Wi - W N , can all be different and are 
typically chosen to be less than the coherence length of the superconducting 
material of terminals 110-1 through 110-N. Terminals 110-1 through 110-N can 


be fabricated from any superconducting material. Two dimensional electron gas 
junction 140 can, for example, be formed of InAs. In an exemplary embodiment 
the two dimensional electron gas is formed of an InAs layer deposited on an AlSb 
substrate. Terminals 110-1 through 1 10-N can be formed of niobium and can be 
deposited on the InAs layer. In accordance with certain embodiments of the 
invention, multi-terminal junction 120 can be used as a circuit element in 
superconducting circuits. 

In general, multi-terminal junction 120 can include any number of 
terminals coupled by any types of junction. The embodiments of junction 120 
shown in Figures 1 A through 1M are exemplary only. One skilled in the art will 
recognize a multitude of variations from those shown above. Those variations are 
intended to be within the scope of this disclosure. 

Additionally, qubit 100 shown in Figure 13 illustratively shows terminals 
1 10-N and 110-1 coupled to portions 124 and 125 of superconducting loop 122. In 
general, superconducting loop 122 can be formed with any arbitrary pair of 
terminals 110-1 through 110-N. Further, qubit 100 shows phase shifter 123 as part 
of superconducting loop 122. 

Figure 2A shows a plan view of an embodiment of a two terminal phase 
shifter 123 having a S/N/D/N/S heterostructure. Phase shifter 123 of Figure 2 A 
includes an s-wave superconducting terminal 210 coupled to a normal metal 
connector 250 which is coupled to a d-wave superconductor 240 which, in turn, is 
coupled to a normal metal connector 25 1 which is coupled to an s-wave 
superconducting terminal 21 1. In some embodiments all lengths and widths L S o, 
L S i, L S 2, L S3 , W S o, and W S i, indicating the termination of terminals 210 and 211, 
can all be different. In some embodiments, the terminal lengths and widths can all 
be less than about five microns. 

Modifying the angle of contact between the d-wave superconductor 240 
and each of the two external terminals 210 and 211 modifies the phase shift 
acquired in transit through phase shifter 123 in a known way. For example, if the 
normal metal connectors, 250 and 251, are at a right angle to each other, the total 
phase is shifted by n through the phase shifter 123. Furthermore, if the normal 


metal connectors were directly opposite (0° apart), then there would be no 
accumulated phase shift. Following from this, any angle between 0° and 90°, 
where the angle is represented by 0, would lead to a phase shift of 20. Figure 2B 
illustrates an exemplary embodiment of a 7u-phase shifter. The angle 0 is 90° and 
5 the normal metal connector 250 is directly parallel with the crystal alignment of the 
junction 240. There is no restriction that the normal metal connector 250 be 
directly parallel with the crystal alignment of the junction 240. 

The physical characteristics of the normal metal connectors 250 and 25 1 
can be chosen so as to provide standard Josephson junction connections between 
1 0 terminal 210 and d-wave superconductor 240 and terminal 211 and d-wave 

superconductor 240, respectively. Currents flowing in terminals 210 and 211 are 
labeled Iso and Isi, respectively. The dimensions of d-wave superconductor 240 
and connectors 250 and 251 are not critical. 

In some embodiments terminals 210 and 211 can be niobium (Nb), 
15 aluminum (Al), lead (Pb) or tin (Sn). In some embodiments d-wave 

superconductor 240 can be YBa 2 Cu 3 07-d, where 0.2<d<0.8. In accordance with an 
exemplary embodiment, terminals 210 and 211 can be made of niobium, 
connectors 250 and 251 of gold, and d-wave superconductor 240 of YBa 2 Cu 3 06.68- 
Lengths Lso, Lsi, Ls2, and Ls3 can be approximately 0.5 microns, widths Wso and 

2 0 Wsi can be approximately 0.5 microns, and connectors 250 and 25 1 can be 

approximately 0.05 microns thick. The embodiment of phase shifter 123 shown in 
Figure 2B will produce a phase shift of n in the superconducting order parameter 
accumulated in transition between terminals 210 and 211. 

Figure 2C shows apian view of another embodiment of two terminal phase 
25 shifter 123. Phase shifter 123 as shown in Figure 2C includes a hetero structure 

containing a grain boundary junction 260 between two lattice-mismatched d-wave 
superconductors 241 and 242. In some embodiments the d-wave superconductors 
241 and 242 can be YBa 2 Cu 3 07-d, where 0.2<d<0.8. Modifying the angle of 
mismatch of the d-wave order parameters at grain boundary 260 between 

3 0 superconductors 241 and 242 affects the phase shift across grain boundary 260 in a 

known way. For example, Figure 2C illustrates a mismatch angle of 45° which 
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would lead to a 7t/2-phse shift. The behavior and resulting phase shift of such 
junctions are well known and are described in detail in C. Binder, A. van Otterlo, 
and G. T. Zimanyi, "Tunnel Junctions of Unconventional Superconductors," Phys. 
Rev. B 51, 12904-07 (1995), and furthermore in R. R. Schultz, B. Chesca, B. 
5 Goetz, C. W. Schneider, A. Schmehl, H. Bielefeldt, H. Hilgenkamp, J. Mannhart, 
and C. C. Tsuei, "Design and Realization of an all d-Wave dc 7t-Superconducting 
Quantum Interference Device," Applied Physics Letters, 76, p. 912-14 (2000), each 
of which are herein incorporated by reference in their entirety. In some 
embodiments the mismatch at grain boundary 260 can be achieved by deposition 

10 of the d-wave superconductor onto a bi-crystal substrate with an existing lattice- 
mismatched grain boundary which the d-wave superconductor inherits. 

Figure 2D shows a cross sectional view of phase shifter 123 as shown in 
Figure 2C. Superconductors 242 and 241 are grown on substrate crystals 271 and 
270, respectively. Substrate crystals 271 and 270 can be mounted on substrate 

1 5 272. In some embodiments bi-crystal substrate 270 and 271 can be an insulator 

such as SrTiC>3 (strontium titanate) or Ti:Al 2 03 (sapphire) which are commercially 
available. 

In some embodiments, grain boundary 260 can be created by using a bi- 
epitaxial method where a d-wave superconductor is deposited onto a substrate 
2 0 containing seed layers upon which the d-wave superconductor grows in a different 
crystallographic direction than on the substrate itself. In some embodiments the 
substrate can be an insulator such as strontium titanate and the seed layers can be 
CeO (cerium oxide) or MgO (magnesium oxide). See F. Tafuri, F. Carillo, F. 
Lombardi, F. Miletto Granozio, F. Ricci, U. Scotti di Uccio, A. Barone, G. Testa, 

2 5 E. Sarnelli, J.R. Kirtley, "Feasibility of Biepitaxial YBa 2 Cu 3 07-x Josephson 

Junctions for Fundamental Studies and Potential Circuit Implementation, Los 
Alamos preprint cond-mat/0010128, published Phys. Rev. B 62, 14431-38 (2000), 
which is herein incorporated by reference in its entirety. 

In some embodiments normal metal connector 250 couples d-wave 

3 0 superconductor 241 to s-wave superconducting terminal 211. In some 

embodiments normal metal connector 251 couples d-wave superconductor 242 to 
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s-wave superconducting terminal 210. In some embodiments normal metal 
connectors 250 and 251 can be gold (Au), silver (Ag), platinum (Pt), or any other 
normal metal substance; and s-wave superconducting terminals 210 and 211 can be 
aluminum (Al), niobium (Nb), or any other superconductor with s-wave pairing 
5 symmetry. 

In some embodiments all lengths and widths Lso> Lsi> Ls2> Ls3> Wso 5 and 
W S i can all be different, in some embodiments each of the lengths can be less than 
about one micron. The physical characteristics of normal metal connectors 250 
and 251 can be chosen so as to provide standard Josephson junction connections 

1 0 between terminals 2 1 0 and d-wave superconductor 241 and terminals 211 and d- 

wave superconductor 240, respectively. Currents flowing in terminals 210 and 211 
are labeled Iso and I S i, respectively. The dimensions of d-wave superconductor 
240 and connectors 250 and 25 1 are not critical. 

In accordance with an exemplary embodiment of phase shifter 123 as 

1 5 shown in Figure 2C, terminals 210 and 211 can be made of niobium, connectors 
250 and 251 of gold, and d-wave superconductor 240 and 241 of YBa 2 Cu 3 06.68. 
Lengths L S o, L S i, L S 2, and L S3 can be approximately 0.5 microns, widths Wso and 
W S i can be approximately 0.5 microns, and connectors 250 and 25 1 can be 
approximately 0.05 microns thick. The d-wave superconductors 240 and 241 can 

2 0 have a symmetric 22.5 / 22.5 degree lattice mismatch, in which the 

crystallographic a-axis of d-wave superconductor 240 makes an angle of +22.5 
degrees with grain boundary 260 and the crystallographic a-axis of d-wave 
superconductor 241 makes an angle of -22.5 degrees with grain boundary 260. 
This type of grain boundary is typically called a symmetric 45 degree grain 

2 5 boundary, as the angle between the crystallographic a-axes of superconductors 240 

and 241 is 45 degrees. This embodiment will produce a phase shift of n in the 
superconducting order parameter accumulated in transition across grain boundary 
260. It is also "quiet" in the sense that no spontaneous supercurrents or magnetic 
fluxes are produced at a symmetric 45 degree grain boundary and therefore noise 

3 0 due to phase shifter 123 in a superconducting electronic circuit is reduced. 
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Figure 2E shows a plan view of another embodiment of a two terminal 
phase shifter 123. Phase shifter 123 of Figure 2E includes a junction between s- 
wave superconductor 210, ferromagnetic region 276 and s-wave superconductor 
211. In this embodiment the s-wave superconductor / ferromagnet / s-wave 
5 superconductor junction is in the axis normal to the plane shown in Figure 2E. 
Figure 2F shows a cross sectional view of phase shifter 123 with ferromagnet 275 
between s-wave superconductor 210 (on the bottom) and s-wave superconductor 
21 1 (on the top). An insulating barrier 275, shown in gray in Figure 2E, provides 
insulation between terminals 210 and 211. 

1 0 Modifying the geometry of the ferromagnetic region 275 can change the 

angle of the phase shift in a known way. In Figure 2E, lengths L S i and L S 3, as 
before, indicate the lengths of terminals 210 and 211, respectively. H TO and H T i 
indicate the distance between the edge of terminals 210 and 21 1, respectively, and 
the edge of insulation region 275. The quantities H F and W F indicate the height 

1 5 and width of ferromagnetic region 276, respectively. The quantity Dti indicates 
the distance between the edge of terminal 21 1 and the edge of terminal 210. In 
some embodiments, lengths and widths D T i, H T i, L S 2, H T o, Wso, and Wsi can be all 
different and, in some embodiments, all be less than about five microns. In some 
embodiments lengths H F and W F can be different and, in some embodiments, can 

2 0 be less than about one micron, with these lengths chosen so as to give a required 
phase shift. Currents flowing in terminals 210 and 21 1 are labeled I S o and I S i, 
respectively. 

In some embodiments terminals 210 and 211 can be niobium (Nb), 
aluminum (Al), lead (Pb), tin (Sn), or any other superconductor with s-wave 

2 5 pairing symmetry. In some embodiments insulating region 275 can be aluminum 

oxide (AIO2) or any other insulating material. In some embodiments 
ferromagnetic region 276 can be an alloy of copper and nickel (Cu:Ni) or any other 
ferromagnetic material. One method for fabricating an example of phase shifter 
123 as shown in Figures 2E and 2F is described in V. V. Ryazanov, V. A. 

3 0 Oboznov, A. Yu. Rusanov, A. V. Veretennikov, A. A. Golubov, J. Aarts, 

"Coupling of Two Superconductors Through a Ferromagnet: Evidence for a n- 
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Junction," Los Alamos preprint server cond-mat/0008364, submitted to Phys. Rev. 
Lett. (2000), which is herein incorporated by reference in its entirety. 

Figure 2G shows a plan view of another embodiment of a two terminal 
phase shifter 123 having a ferromagnetic region 276 imbedded in a junction 
5 between s-wave superconductors 210 and 2 1 1 . In this embodiment the s-wave 
superconductor / ferromagnet / s-wave superconductor junction is in the plane of 
Figure 2G. Thus, the ferromagnetic region, 276, is directly in the plane with of the 
terminals 210 and 21 1. Modifying the geometry of the ferromagnetic region 276 
can change the angle of the phase shift in a known way. In some embodiments 

1 0 lengths and widths D T i, H T i, L S2 , W S o, and W S i can be all different and, in some 
embodiments, all be less than about five microns. In some embodiments lengths 
H F and W F can be different and less than about one micron, with these lengths 
chosen so as to give a required phase shift. Currents flowing in terminals 210 and 
21 1 are labeled I S o and I S i, respectively. In some embodiments terminals 210 and 

15 2 1 1 can be niobium (Nb), aluminum (Al), lead (Pb) tin (Sn), or any other 
superconductor with s-wave pairing symmetry. In some embodiments 
ferromagnetic region 276 can be an alloy of copper and nickel (Cu:Ni) or any other 
ferromagnetic material. Ferromagnetic region 276 can be prepared by, for 
example, implantation of a ferromagnetic substance into a superconducting 

2 0 junction. 

Figure 3A shows a plan view of an embodiment of a four-terminal qubit 
100 having junction 120 and phase shifter 123 in superconducting loop 122. 
Junction 120 can be any four-terminal junction. Examples of embodiments of 
junction 120 are shown in Figures 1 A through 1M. Examples of embodiments of 

2 5 phase shifter 123 are shown in Figures 2 A through 2G. In Qubit 1 00 as shown in 

Figure 3 A, terminal 110-1 can be reduced from width Wi to width Wji before 
being coupled to junction 120, terminal 110-2 can be reduced from width W 2 to 
width W J2 before being coupled to junction 120, terminal 110-3 can be reduced 
from width W 3 to width W J3 before being coupled to junction 120, and terminal 

3 0 110-4 can be reduced from width W 4 to width W J4 before being coupled to junction 

120. Terminal 210 of phase shift 123 is coupled to portion 124 having width W S 4 
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and terminal 21 1 of phase shift 123 is coupled to portion 125 having width W S6 - 
Terminal 110-3 is coupled to portion 125 having width W S 5- 

Terminals 110-1, 110-2, 110-3, 110-4, 210, and 211 along with portions 
124 and 125 can be made of any superconducting material compatible with the 
particular embodiments of phase shifter 123 and four-terminal junction 120. 
Superconducting currents Ii, h, h and Lt can exist in terminals 1 10-1, 1 10-2, 1 10-3 
and 110-4, respectively. Widths Wji through W J4 are constrained only by the 
requirement of compatibility with junction 120 and widths W S o and W S i are 
compatible with phase shifter 123. In some embodiments, widths Wi through W 4 , 
Wji through W J4 , and W S o through W S 6 can all be less than about 10 microns. 
Lengths Li, L 2 , Dji, D J2 , D J3 , D J4 , D P0 , and D P i are all compatible with phase 
shifter 123 and four-terminal junction 120 and, in some embodiments, can be all 
less than about 10 microns. Superconducting loop 122 can be threaded by a 
magnetic flux O, which may contain contributions from a spontaneous 
supercurrent in the loop and externally applied magnetic fields. 

Four-terminal junction 120 couples one end of each of the four terminals 
110-1 through 110-4 by, for example, constriction junctions, tunnel junctions, two- 
dimensional electron gas structures, or combinations of these. The choice of the 
physical sizes of the elements in four-terminal junction 120 that couple the four 
terminals 110-1 through 110-4 also affects the function of qubit 100. To achieve a 
small magnetic flux O in superconducting loop 122 (which is desirable for 
coherence consideration) and maximum influence of the transport current It=Ii=I2 
on the properties of superconducting loop 122, in some embodiments the links in 
the transport loop are much wider than the ones in superconducting loop 122, that 
is, Wji and W J2 are much larger than W J3 and W J4 . A small magnetic flux O can 
exist even in the absence of external magnetic fields because of spontaneous 
supercurrents. In those embodiments, the height of the potential energy barrier 
between the two degenerate quantum states of the qubit quantum system of qubit 
100 will be affected most pronouncedly by the applied transport current It=Ii=I 2 . 

SQUID loop 122 with intrinsic phase shifter 123 can provide a basic block 
for construction of qubit 100 but can also be utilized for demonstration of 
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macroscopic quantum tunneling and incoherent quantum noise in a solid state 
system. As described further below, the macroscopic quantum tunneling in a set of 
independent four-terminal qubits 100 with intrinsic phase shifters 123 (i.e., with no 
entanglements between individual qubits) permits construction of a random 
5 number generator that generates random series with zero correlation between 
numbers in the random series. 

Four-terminal qubit 100 with intrinsic phase shifter 123 includes a 
superconducting loop 122 linking two of terminals 110-1 through 110-4, terminals 
1 10-4 and 110-3 in Figure 3 A. Four-terminal qubit 100 with intrinsic phase shifter 

10 123 can be made superconducting by reducing the temperature of qubit 1 00 below 
the superconducting critical temperature T c of all of the superconducting materials 
utilized in the formation of qubit 100. Four-terminal junction 120 maybe either 
symmetric or asymmetric, as was discussed with respect to Figures 1 A through 
1M. The superconducting materials from which qubit 100 is constructed are 

1 5 constrained only by the requirement of compatibility with phase shifter 123 and 
four-terminal junction 120 and otherwise may have any pairing symmetry. For 
example, materials used may be s-wave, for example, niobium or aluminum, or d- 
wave, such as a high-Tc cuprate, for example YBa 2 Cu307- x , or any 
superconducting material in which the Cooper pairs are in a state with non-zero 

2 0 orbital angular momentum. 

Four-terminal qubit 100 with intrinsic phase shifter 123 can be formed on 
an insulating substrate such as, for example, strontium titanate or sapphire. (See 
Figure 2C, for example). Phase shifter 123 may be any structure that shifts the 
phase of the superconducting order parameter in transition across the structure. 

2 5 Examples of phase shifter 123 are shown in Figures 2 A through 2G. 

Superconducting loop 122 can include multi-crystalline d-wave superconducting 
material where the phase shift is caused by transition through the crystalline 
boundaries of the multi-crystalline material. In such embodiments, phase shifter 
123 is distributed throughout portions 124 and 125 of Figure 3A. 

3 0 Figure 3B shows a plan view of an embodiment of a N-terminal qubit 100 

according to the present invention. Qubit 100 includes terminals 110-1 through 
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1 10-N where superconducting loop 122 is formed between arbitrary terminals 1 10- 
I and 110-K. Terminals 1 10-1 and 1 10-K can be any pair of terminals 110-1 
through 1 10-N. Qubit 100 includes phase shifter 123 in superconducting loop 122 
between portion 125 and portion 124. Multi-terminal junction 120 couples 
terminals 110-1 through 1 10-N. Examples of embodiments of phase shifter 123 
are shown in Figures 2 A through 2G. Examples of embodiments of junction 120 
are shown in Figures 1A through 1M. Terminals 110-1 through 1 10-N can be 
made of any superconducting material compatible with the particular embodiment 
of phase shifter 123 and N-terminal junction 120. Terminals 110-1 through 1 10-N 
carrying superconducting currents Ii through I N , respectively. The transport 
current I T is the current directed to terminal 1 10-1, Ii. 

Widths Wi through W N> widths Wji through Wjn, and widths W S o through 
W S 6 are only constrained by the requirement of compatibility with phase shifter 
123 and N-terminal junction 124 and, in some embodiments, typically all less than 
about 10 microns. Terminals 110-1 through 1 10-N of qubit 100 as shown in 
Figure 3B include a narrowing of the width of the terminal as the terminal is 
couple to junction 120. Additionally, portions 124 and 125 are narrowed to 
coupled with terminals 210 and 21 1 of phase shifter 123. Lengths Li through L N , 
Dji through Djm, L so , and L S i are only constrained by the requirement of 
compatibility with phase shifter 123 and N-terminal junction 120 and, in some 
embodiments, are typically all less than about 10 microns. The superconducting 
loop can be threaded by a magnetic flux O, which may contain contributions from 
a spontaneous supercurrent in the loop and externally applied magnetic fields. 

N-terminal junction 120 couples one end of each of terminals 110-1 
through 1 10-N by, for example, constriction junctions, tunnel junctions, two- 
dimensional electron gas structures, or combinations of these. The choice of the 
physical sizes of the elements in N-terminal junction 120 also affects the function 
of qubit 100. To achieve a small magnetic flux O in superconducting loop 122 
(which is desirable for coherence consideration) and maximum influence of the 
transport currents Ii through I N on the properties of superconducting loop 122, in 
some embodiments the terminals in the transport terminals (i.e., terminals 110-1 


through 1 10-N that are not terminals 1 10-1 and 1 10-K) are much wider than 
terminals 1 10-1 and 1 10-K at junction 120, that is, Wji through Wjn, excluding Wji 
and W JK are much larger than Wji and W JK - A small magnetic flux O can exist, 
even in the absence of external magnetic fields, because of spontaneous 
5 supercurrents. In those embodiments, the height of the potential energy barrier 
between the two degenerate quantum states of the qubit quantum system will be 
affected most pronouncedly by the applied transport currents Ii through I N . 

N-terminal qubit 100 with intrinsic phase shifter 123 can provide a basic 
block for construction of a qubit but can also be utilized for demonstration of 

1 0 macroscopic quantum tunneling and incoherent quantum noise in a solid state 

system. As described further below, the macroscopic quantum tunneling in a set of 
independent N-terminal qubits with intrinsic phase shifters (i.e., with no 
entanglements between individual qubits) permits construction of a random 
number generator that generates random series with zero correlation between 

1 5 numbers in the random series. 

N-terminal qubit 100 with intrinsic phase shifter 123 includes a 
superconducting loop 122 linking two of the N terminals, terminals 1 10-1 and 110- 
K. N-terminal qubit 100 with intrinsic phase shifter 123 is made superconducting 
by reducing the temperature of qubit 100 below the superconducting critical 

2 0 temperature T c of all of the superconducting materials in qubit 100. N-terminal 
junction 120 may be either symmetric or asymmetric. The superconducting 
materials from which qubit 100 A is constructed are constrained only by the 
requirement of compatibility with phase shifter 123 and N-terminal junction 120 
and otherwise may have any pairing symmetry. For example, materials used may 

2 5 be s-wave, for example, niobium or aluminum, or d-wave, such as a bigh-Tc 

cuprate, for example YBa2Cu 3 07- x , or any superconducting material in which the 
Cooper pairs are in a state with non-zero orbital angular momentum. 

N-terminal qubit 100 with intrinsic phase shifter 123 can be formed on an 
insulating substrate such as, for example, strontium titanate or sapphire. Phase 

3 0 shifter 123 may be any structure that shifts the phase of the superconducting order 

parameter in transition across the structure. Examples of embodiments of phase 
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shifter 123 are shown in Figures 1A through 1G. Additionally, phase shifter 123 
can be incorporated into portions 124 and 125 if portions 124 and 125 are of a 
multi-crystalline d-wave superconducting material, where the phase shift is caused 
by transition through the crystalline boundaries of the multi-crystalline material. 

Figure 3C shows a schematic diagram of qubit 100. Terminals 1 10-1 
through 1 10-N are coupled at one end through junctions 140-1 through 140-(N-1) 
in multi-terminal junction 120. Terminals 110-1 and 1 10-K are coupled to portions 
124 and 125. Phase shifter 123 is included in superconducting loop 122. 

Figure 4A shows a plan view of an embodiment of a qubit array 400. Qubit 
array 400 includes a series of multi-terminal qubits 401-1 through 401-M, each of 
which includes an intrinsic phase shifter 402-1 through 402-N coupled in series to 
terminals 410 and 411, which are themselves coupled to form a loop. As such, 
qubit array 400 as shown in Figure 4A is an example of a qubit register structure. 
Multi-terminal qubits 400-1 through 400-N can be embodiments of qubit 100 as 
shown in Figures 3A through 3C. Each of qubits 401-1 through 401-2 can be 
threaded by a magnetic flux O 480 -i through O 480 -n, respectively. The loop formed 
by terminals 410 and 422 can be threaded by magnetic flux 0 48 i. Terminals 410 
and 41 1 can be made of any superconducting materials that can be coupled to the 
terminals of qubits 401-1 through 401-N. 

Terminal 410 is coupled to one of the terminals of qubit 401-1. Another of 
the terminals of qubit 401-1 is coupled to a terminal of qubit 401-2. Each of qubits 
401-2 through 401-N are coupled to a terminal of qubit 401-1 through 401-(N-1). 
Additionally, a second terminal of qubit 401-N is coupled to terminal 41 1. 
Connectors joining qubits 401-1 through 401-N can be made of any 
superconducting material compatible with qubits 401-1 through 401-N. 
Superconducting current I L can exist in the loop formed by joining terminals 410 
and 41 1 . Widths W L i through Wi^sr+i) and W SL i through W SL 3 are constrained only 
by the requirement of compatibility with qubits 400-1 through 400-N and are 
typically all less than about 10 microns. Lengths L S li through L S u and L L] 
through Ll^+i) are not critical but, in some embodiments, are typically all less than 
about 10 microns. 


The magnetic flux ® 48 i, which threads the loop formed by joining terminals 
410 and 41 1, may contain contributions from a spontaneous supercurrent in the 
loop and contributions from externally applied magnetic fields. The magnetic 
fluxes O 4 8o-i through O 48 o-N, which thread the loops in multi- terminal qubits 401-1 
5 through 401-N, respectively, may contain contributions from a spontaneous 

supercurrent in the loop and contributions from externally applied magnetic fields. 

In accordance with embodiments of aspects of the present invention, qubit 
array 400 as shown in Fig. 4A demonstrates coupling of different qubit quantum 
systems (i.e., qubits 401-1 through 401-N) and reading out information about this 

1 0 coupling. In embodiments where there is no externally applied magnetic field, the 
magnetic flux <3> 4 8i threading the loop formed by joining terminals 410 and 411 is a 
known function of the quantum states of the qubits 401-1 through 401-N. A 
discussion of the relationship between the magnetic flux 0 4 gi and the quantum 
states of qubits 401-1 through 401-N is given in Appendix A, herein incorporated 

15 by reference in its entirety. 

Therefore, measurement of the flux 0 48 i, for example by an external 
measuring instrument such as a magnetic force microscope, scanning SQUID 
microscope or scanning Hall probe, provides information on the quantum states of 
qubits 401-1 through 401-N. In addition, application of time-dependent external 

2 0 magnetic fields to the superconducting loops of SQUID qubits 401-1 through 401- 
N and to the loop formed by joining terminals 410 and 41 1 can function as an 
operating system to perform specific algorithms using the qubit register 
architecture of qubit array 400 as disclosed in Fig.4A. The relation of the time- 
dependent magnetic field application and the algorithm performed is described in 

2 5 Appendix A. 

Figure 4B shows a plan view of an embodiment of a qubit array 400. Qubit 
array 400 includes series coupled qubits 401-1 through 401-N coupled to terminals 
410 and 411. Each of qubits 401-1 through 401-N includes a superconducting loop 
122 with intrinsic phase shifter 402-1 through 402-N, respectively. Terminals 410 

3 0 and 411 are coupled to an external source of transport current to provide current It- 

Multi-terminal qubits 401-1 through 401-N are embodiments of aspects of the 


invention described in Figs.3A through 3C and can be any multi-terminal qubit. 
Qubits 401-1 through 401-N can be threaded by magnetic fluxes <I>480-i through 
®480-n, respectively. 

Terminals 41 0 and 41 1 can be made of any superconducting material 
5 compatible with the choice of qubits 40 1 - 1 and 40 1 -N. Connectors j oining qubits 
401-1 through 401-N can be made of any superconducting material compatible 
with the choice of qubits 401-1 through 401-N. Superconducting current I T , which 
can be present in terminals 410 and 411, typically arises from an external transport 
current source. 

1 0 Widths W L i through W L (n+i> and W S li through W S L3 are only constrained 

by the requirement of compatibility with qubits 401-1 through 401-N and, in some 
embodiments, are typically all less than about 10 microns. Lengths L L i through 
Ll(n+i>, Lsli and L S l3 are not critical but are typically all less than about 10 
microns. The magnetic fluxes <E>4so-i through O480-N which thread the 

15 superconducting loops in qubits 401-1 through 401-N, respectively, may contain 
contributions from a spontaneous supercurrent in the loop and contributions from 
externally applied magnetic fields. 

In accordance with an embodiment of an aspect of the invention, qubit 
array 400 of Fig.4B demonstrates an apparatus for coupling different qubit 

2 0 quantum systems and reading out information about this coupling. In the case 
where there is no externally applied transport current (i.e., no transport current 
through other ones of the terminals of qubits 401-1 through 401-2 other than those 
coupled to terminals 410 and 41 1), the current It in terminals 410 and 41 1 is a 
known function of the quantum states of qubits 401-1 through 401-N, as discussed 

2 5 in Appendix A. Therefore measurement of the current I T , for example by an 

external measuring instrument such as a single electron transistor, provides 
information on the logical states of qubits 401-1 through 401-N. In addition, 
application of time-dependent externally applied transport current Ix(t) to terminals 
410 and 411 and time-dependent external magnetic fields Bi(t) through B N (t) to the 

3 0 superconducting loops of qubits 40 1 - 1 through 40 1 -N, respectively, can function 

as an operating system to perform specific algorithms using the qubit register 
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architecture of qubit array 400 shown in Fig.4B. The relation of the time- 
dependent magnetic field B(t) and transport current I T (t) application and the 
algorithm performed is described in Appendix A. 

Figure 5 shows apian view of an qubits 500 and 501 coupled at junction 
120 where junction 120 is, in this example, a six-terminal junction. Qubits 500 and 
501, then, share one junction 120. Qubit 500 includes a superconducting loop 
coupled into junction 120 with terminals 110-5 and 1 10-6. Qubit 501 includes a 
superconducting loop coupled into junction 120 with terminals 1 10-2 and 110-3. 
Terminals 110-1 and 110-4 can be coupled to a source of transport current I T . 
Further, qubits 500 and 501 include phase shifters 523-1 and 523-2, respectively. 

Examples of embodiments of phase shifters 523-1 and 523-2 are shown as 
phase shifter 123 in Figures 2A through 2G. Examples of embodiments of 
junction 120 are shown in Figures 1A through 1M. Examples of embodiments of 
qubits 500 and 501 are shown in Figures 3a through 3C. 

The superconducting loops of qubits 500 and 501 can be threaded by 
magnetic fluxes O 580 and <& 5 8i, respectively. Terminals 110-1 through 1 10-6 can 
be made of any superconducting material compatible with the choice of qubits 500 
and 501. Superconducting current I T in terminals 110-1 and 110-4 typically arises 
from an external transport current source. The superconducting loops of qubits 
500 and 501 carry superconducting currents I Q0 and Iqj, respectively. The physical 
dimensions of terminals 110-1 through 1 10-6 are constrained only by the 
requirement of compatibility with qubits 500 and 501 and, in some embodiments, 
are typically all less than about 10 microns. The magnetic fluxes 3> 58 oand <D 58 i 
which thread the superconducting loops of qubits 500 and 501, respectively, may 
contain contributions from a spontaneous supercurrent in the loop and 
contributions from externally applied magnetic fields. 

In accordance with an embodiment of an aspect of the invention, qubits 500 
and 501 as shown in Fig. 5 demonstrates an apparatus for coupling two different 
qubit quantum systems (i.e., the quantum systems of qubits 500 and 501) and 
reading out information about this coupling. In the case where there is no 
externally applied transport current, the current I T in terminals 110-1 and 1 10-4 is a 


known function of the quantum states of qubits 510 through 51 1, as discussed in 
Appendix A. Therefore measurement of the current I T , for example by an external 
measuring instrument such as a single electron transistor, provides information on 
the quantum states of qubits 500 and 501. In addition, application of time- 
dependent externally applied transport current I T (t) to terminals 110-1 and 1 10-4 
and time-dependent external magnetic fields B(t) to the superconducting loops of 
qubits 500 through 501 can function as an operating system to perform specific 
algorithms using the qubit register architecture qubits 500 and 501 as shown in 
Fig. 5. The relation between the time-dependent magnetic field and transport 
current application and the algorithm performed is described in Appendix A. 

Figure 6A shows a plan view of an embodiment of a qubit array 610. Qubit 
array 690 includes qubits 600-1 through 600-(N+l) coupled by six-terminal 
junctions 601-1 through 601-N. Examples of embodiments of junctions 601-1 
through 601-N are shown as junction 120 in Figures 1A through 1M. The 
superconducting loop of each of qubits 600-2 through 600-N includes two of 
junctions 600-1 through 600-N whereas the superconducting loops of qubits 600-1 
and 600-(N+l) includes one of junctions 600-1 and 600-N. Each of Junctions 601- 
1 through 601-N, as shown in Figure 6A, includes six terminals labeled 1 through 
6 in counter-clockwise fashion. In an arbitrary one of junctions 601-1 through 
601-N, terminals 1 and 4 are coupled to receive transport currents I T1 through I™, 
respectively. Terminals 2 and 3 form part of a superconducting loop (e.g., 
terminals 2 and 3 of junction 601-2 forms part of the superconducting loop of qubit 
600-3). Terminals 5 and 6 form part of a separate superconducting loop (e.g., 
terminals 5 and 6 of junction 601-2 forms part of the superconducting loop of qubit 
600-2). The superconducting loop of qubits 600-1 through 600-(N+l) each include 
a phase shifter 602-1 through 602-(N+l), respectively. Examples of embodiments 
of phase shifter 602-1 through 602-(N+l) are shown as phase shifter 123 in Figures 
2A through 2G. 

The superconducting loop of qubits 600-1 through 600-(N+l) can be 
threaded by magnetic fluxes <D 6 80-i through O 68 o-(N+1), respectively. The 
superconducting loops of qubits 600-1 through 600-N can include any 


superconducting material compatible with junctions 601-1 through 601-N and 
phase shifters 602-1 through 602-(N+l). Terminals 1 and 4 of each of junctions 
601-1 through 601-N can carry superconducting currents I T i through Itn, which 
typically arise from external transport current sources. The superconducting loops 
5 of each of qubits 600-1 through 600-(N+l) can carry superconducting currents 
Lsl-i through L S l-(n+i> respectively. All length and width scales are constrained 
only by the requirement of compatibility with junctions 601-1 through 601-N and 
phase shifters 602-1 through 602-(N+l), as has been previously discussed. The 
magnetic fluxes d> 6 8o-i through O 680 -(n+i) which thread the loops in multi-terminal 
1 0 qubits 600 through 600-(N+l) respectively may contain contributions from a 
spontaneous supercurrent in the loop and contributions from externally applied 
magnetic fields. 

In accordance with an embodiment of an aspect of the invention, qubit 
array 610 as shown in Fig.6A demonstrates an apparatus for coupling a plurality of 

15 different qubit quantum systems and reading out information about the quantum 
states of the qubit quantum systems. In the case where there are no externally 
applied transport currents, the currents Iti through Itn in terminals carried by 
terminals 1 and 4 of junctions 601-1 through 601-N are known functions of the 
logical states of the qubits 600-1 through 600-(N+l), as described in Appendix A. 

2 0 Therefore measurement of the currents In through Itn, for example by external 

measuring instruments such as single electron transistors, provides information on 
the logical states of the qubits 600-1 through 600-(N+l). In addition, application 
of time-dependent externally applied transport currents I T i(t) through IinOO to 
terminals 1 and 4 of junctions 601-1 through 601-N, respectively, and time- 

2 5 dependent external magnetic fields to the superconducting loops of qubits 600-1 

through 600-(N+l) can function as an operating system to perform specific 
algorithms using the qubit register architecture of qubit array 610 as shown in 
Fig. 6 A. The relation of the time-dependent magnetic field and transport current 
application and the algorithm performed is described in Appendix A. Readout 605 

3 0 reads the quantum states of each of qubits 600-1 through 600-(N+l). 
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Figure 6B shows a plan view of an embodiment of a qubit register 680 
from an array of M -qubit registers 690-1 through 690-M. An example of an 
embodiment of qubit registers 690-1 through 690-M is shown as qubit array 610 of 
Figure 6 A. Registers 690-1 through 690-M are coupled together through junctions 
5 695-1 1 through 695-(M-l)N. With reference to Figure 6A, terminals 1 and 4 of 
each of junctions 601-1 through 601-N for a particular one of registers 690-1 
through 690-M are coupled to counterpart terminals 4 and 1, respectively, of 
junctions 601-1 through 601-N adjoining ones of registers 690-1 through 690-M. 
In the embodiment shown in Figure 6B, each of junctions 695-1,1 through 695-(M- 

10 1 ),N are four-terminal junctions, examples of which are given as junction 1 20 in 
Figures 1A through 1M. In Figure 6B, the four terminals of each of 695-1,1 
through 695-(M-l),N are labeled counterclockwise as terminals 1-4, as shown for 
junction 695-1,1. As an example, terminal 3 of junction 695-1,1 is coupled to 
terminal 1 of junction 601-1 of array 690-1 and terminal 1 of junction 695-1,1 is 

1 5 coupled to terminal 4 of junction 60 1 -1 of array 690-1 . Terminal 2 is coupled to 
terminal 4 of junction 695-1,2. Terminal 4 of junctions 695-1,1 through 695-(M- 
1),1 can be coupled to current sources to receive currents Iai through Ia(m-i), 
respectively. Further, terminal 1 of each of junctions 601-1 through 601-N of array 
690-M can be coupled to receive currents Ibi through Ibn 5 respectively. All length 

2 0 and width scales are constrained only by the requirement of compatibility with 
qubit registers 690-1 through 690-M as described in Fig. 6 A. 

Magnetic fluxes ® 6 8o-u through 0 6 8o-(M),(N) can be embraced by 
superconducting loops contained in qubit registers 690-1 through 690-M. The 
magnetic fluxes <I>680-ii through <I>68o-(mxn)b which thread the superconducting loops 

2 5 in qubits 600-1 through 600-N of each of qubit registers 690-1 through 690-M, 

respectively, may contain contributions from a spontaneous supercurrent in the 
loop and contributions from externally applied magnetic fields. 

In accordance with an embodiment of an aspect of the invention, the array 
of registers 690-1 through 690-M shown in Fig.6B demonstrates an apparatus for 

3 0 coupling a plurality of different qubit quantum systems and reading out 

information about the logical states of the qubit quantum systems. In the case 
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where there are no externally applied transport currents, the currents I A i through 
I A (M-i) and currents I B i through I B (N-i) are known functions of the quantum states of 
the qubits in qubit registers 690-1 through 690-M, as discussed in Appendix A. 
Therefore measurement of the currents I A i through Ia(m-i) and I B i through Ib(n-i>, 
5 for example by external measuring instruments such as single electron transistors, 
provides information on the quantum states of the qubit 600-1 through 600-N in 
each of qubit registers 690-1 through 690-M. In addition, application of time- 
dependent externally applied transport currents I A i(t) through lA(M-i)(t) to terminal 
4 of junctions 695-11 through 695-(M-l)l, respectively, and iBi(t) through lB(N-i)(t) 

10 to terminal 1 of junctions 601-1 through 601-N of qubit register 690-M, 
respectively, and time-dependent external magnetic fields B(t) to the 
superconducting loops of qubits 600-1 through 600-N of each of qubit register 
690-1 through 690-M function as an operating system to perform specific 
algorithms on qubit register 680 as shown in Fig.6B. The relation of the time- 

15 dependent magnetic field and transport current application and the algorithm 
performed is described in Appendix A. 

In accordance with another embodiment of an aspect of the present 
invention, Fig.7 demonstrates an apparatus for measuring voltages 700. Apparatus 
700 is a radio-frequency single electron transistor electrometer and is well-known 

2 0 and described, for example, in A. N. Korotkov and M. A. Paalanen, "Charge 

Sensitivity of Radio-Frequency Single Electron Transistor, Appl. Phys. Lett. 74, 26 
(1999), which is herein incorporated by reference in its entirety. Apparatus 700 
can be utilized in readout 605 of Figure 6A to read the quantum state of each of 
qubits 600-1 through 600-(N+l). 

2 5 The single-electron transistor (SET) 709 can be made of any material that 

displays a Coulomb blockade effect, for example niobium, aluminum, lead, tin, 
and any high-temperature superconducting cuprate. P. Joyez, P. Lafarge, A. Filipe, 
D. Esteve, and M. H. Devoret, "Observation of Parity-Induced Suppression of 
Josephson Tunneling in the Superconducting Single Electron Transistor", Phys. 

3 0 Rev. Letters, Vol. 72, No. 15, 2458-61 (April 11, 1994), describes operation and 

manufacture of single electron transistors and is incorporated by reference herein 
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in its entirety. SET 709 is placed in a high quality factor tank circuit 712 tuned to 
resonance. Tank circuit 712 includes inductor 707 and capacitor 708. Capacitor 
708 is coupled across SET 709. A third terminal of SET 709 is coupled to 
electrode 710. A radio-frequency or microwave signal 704 is introduced into the 
5 circuit 712. The reflected signal 705 is a strong function of the voltage difference 
between electrode 710 and ground 711. Analysis of reflected signal 705 using 
established techniques allows measurement of the voltage difference between 
electrode 710 and ground 711. 

Read-out of the state of the qubit quantum system may be done via the use 

10 of a single electron transistor (SET) 709 according to known procedures, for 
example, as described in R. J. Schoelkopf, P. Wahlgren, A. A. Kozhevnikov, P. 
Delsing, and D. E. Prober, "The Radio-Frequency Single-Electron Transistor (RF- 
SET): A Fast and Ultrasensitive Electrometer," Science, Vol. 280, 1238-42 (May 
22, 1998). SET 709 may be coupled to three devices (e.g., terminals 710, 71 1 and 

15 712). An electron or Cooper pair can tunnel onto SET 709 when SET 709 is 
uncharged. However, SET 709 is small enough that once an electron or Cooper 
pair tunnels onto SET 709, the charging of SET 709 electrically repels and 
prevents further tunneling onto SET 709. A terminal 710 associated with SET 709 
can change the voltage of SET 709 and de-tune tank circuit 712, changing the 

2 0 characteristics of the reflected wave 705. 

In operation, in order to measure a current, for example one of currents I A i 
through Ia(m-i) or Ibi through Ib(n-i> shown in Figure 6B, electrode 710 is coupled 
to terminal 4 of junctions 495-1 1 through 495-(M-l)N or terminal 1 of junctions 
601-1 through 601-N of qubit register 690-M. The current at those terminals can 
25 be measured by applying signal 704 and monitoring signal 705 of Figure 7. 

All structures in Figures 1 A through 7 can be formed using conventional 
fabrication techniques. Elemental s-wave superconductors may be purchased 
from, for example, CIL Cambridge Isotopes Laboratories Inc. and substrates upon 
which superconducting structures can be formed may be purchased, for example, 

3 0 from Kagaku Gijutsu-sha of Tokyo, Japan or from Shinkosha, Ltd. (c/o Nikko 

Trading Co.). The fabrication process for a four-terminal structure with s-wave 


superconductors is developed in detail, for example, in B. J. Vleeming, "The Four- 
Terminal SQUID," Ph.D. Dissertation, Leiden University, The Netherlands, 1998 
and references therein, which has previously been incorporated into this disclosure 
by reference. 

5 Some embodiments of the invention allow all of the operations that are 

required for quantum computing to be done without the application of external 
magnetic fields. Operations such as read and initialization, as well as operating 
system gates such as application of Pauli operators ct x and cr z , and furthermore 
operations for maintaining coherence in the state of the qubit can be achieved, 

1 0 without the use of external magnetic fields. 

In some embodiments, as illustrated in Figure 8, a qubit 100 comprises a 
five terminal junction 120 and a superconducting loop 122 with phase shifter 123. 
The embodiment of qubit 100 shown in Figure 8 includes terminals 110-1 through 
110-5, where terminals 110-4 and 110-5 are connected to form superconducting 

1 5 loop 122 by a phase shifter 123. Terminals 110-1 through 1 10-3 can be utilized for 
the application of transport current as well as for voltage measurements in reading 
the quantum state of qubit 100. Five terminal junction 120 can be any junction for 
coupling terminals 110-1 through 110-5, such as those discussed with Figures 1A 
through 1M. For example, junction 120 can be a two dimensional electron gas 

2 0 junction and terminals 1 10-1 through 1 10-5 can be of niobium. A controller 800 
can be coupled with terminals 110-1 through 1 10-3 in order to supply currents to 
and measure voltage across various ones of terminals 1 10-1 through 110-3, as 
discussed below. Qubit 100 as shown in Figure 8, then, can be both symmetric and 
asymmetric. 

2 5 Although qubit 100 in Figure 8 is being discussed as a five-terminal qubit, 

one skilled in the art will recognize that methods and examples given here are 
extendable to qubit 100 having any number of terminals 110-1 through 1 10-N. 

The quantum state of the quantum system of qubit 100 of Figure 8, for 
example, can be initialized by controller 800 passing a transport current I T through 

3 0 junction 120 for a sufficient duration, as is further discussed in Appendix A. The 

transport current I T through junction 120 has the effect of biasing the energy states 
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of the quantum system of qubit 100 towards a particular energy state and of 
breaking the ground state degeneracy. Given a sufficient period of time, the 
quantum system of qubit 100 will relax into the lower energy state. By 
manipulating the direction of the transport current I T passing through junction 120, 
qubit 100 can be initialized such that the quantum state of qubit 100 is the desired 
state. For example, controller 800 can apply a bias current from terminal 1 10-1 to 
terminal 110-3, thereby selects a particular state of the quantum system of qubit 
100. By reversing the current through the same terminals, controller 800 can select 
the opposite state of the quantum system. 

A read operation on qubit 100 of Figure 8 can be accomplished based on 
the fact that each of the two degenerate ground states of the quantum system of 
qubit 100 exhibits a unique current-voltage curve with respect to current flowing 
between terminals 110-1 and 110-3. Each of the two degenerate ground states 
results in a different critical current in junction 120. The critical current is the 
current which, if exceeded, results injunction 120 developing a resistance. 
Therefore, determining which of the two critical currents is appropriate for 
junction 120 differentiates between the two degenerate ground states of the 
quantum system. 

The quantum state of the quantum system of qubit 100 of Figure 8, for 
example, can be read by controller 800 passing a transport current I T through 
junction 120 (for example between terminals 110-1 and 110-3). The critical 
current I c of junction 120 is dependent on the quantum state of the quantum system 
of qubit 100, with one state corresponding to a lower value of critical current in 
junction 120 and the opposite state corresponding to a higher value of critical 
current. In any Josephson junction, if the critical current I c is exceeded, dynamical 
effects result and a resistance becomes present in the junction. Thus, determining 
the state of the quantum system of qubit 100 can be accomplished by discerning 
the value of the critical current Ic in the junction, see Appendix A. 

In one method of measuring the quantum state of the quantum system of 
qubit 100, controller 800 applies a transport current It to junction 120 which is 
between the known upper and lower critical current values (i.e., between the values 


of the critical current Ic for each of the quantum states). The upper and lower 
values of the critical current Ic is dependent upon the particular embodiment of 
qubit 100. When the transport current I T is applied, if the system occupies the state 
associated with the lower critical current, then the transport current will have 
5 exceeded the critical current value of the junction, thus resulting in a junction 

resistance, and a corresponding voltage across the terminals (for example, between 
terminals 110-1 and 110-3 to which the transport current I T is applied). 
Alternatively, if the system occupies the high critical current state, no voltage 
across the terminals will result. Controller 800, then, can determine the quantum 

1 0 state of the quantum system of qubit 1 00 by monitoring the voltage across junction 
120 while applying the transport current I T through junction 120. For example, by 
applying a transport current from terminal 110-1 through terminal 1 10-3 of qubit 
100 of Figure 8, and measuring the voltage between terminals 1 10-1 and 110-3, it 
is possible for controller 800 to determine the state of the system, where the 

1 5 presence of a voltage indicates one of the states and the absence of voltage 

represents the opposite state. In some embodiments, controller 800 blocks the flow 
of current through terminal 110-2 by, for example, shorting terminal 110-2. 

A phase gate operation ct z can be performed, for example, on qubit 1 00 of 
Figure 8 by controller 800 applying a transport current It pulse through junction 

2 0 120. In essence, a quantum gate operation operates to bias a particular state of the 
quantum system of qubit 100. Thus, in action, the current pulse producing the a z 
operation can be similar to the read operation discussed above, except that the 
magnitude of the current pulse applied is much less than the magnitude used for 
the exemplary read operation so that no dynamical effect results from its 

2 5 application, see Appendix A. The magnitude and duration of the transport current 

It pulse necessary to affect the <j z operation is specific to the particular 
embodiment of qubit 100. As long as the magnitude of the pulse of transport 
current It is small, application of the transport current It by controller 800 will not 
destroy the quantum superposition of states in the quantum system of qubit 100, 

3 0 but merely weight one of the states as desired. 
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A phase gate operation cr x can be performed, for example, on qubit 100 of 
Figure 8 by application of a transport current I T to terminal 1 10-2 and allowing it 
to escape through both of terminals 110-1 and 110-3. In performing the ct x 
operation, then, controller 800 applies a transport current to junction 120 through 
terminal 1 10-2, and having current flow out of junction 120 through terminals 110- 
1 and 110-3. Application of the transport current in this manner results in 
manipulating the height of the potential barrier that separates the two degenerate 
states of the quantum system of qubit 100, see Appendix A. The magnitude of the 
transport current applied is dependent on the actual configuration of qubit 100. 
When applied for a short duration of time, the height of the potential barrier 
between the two degenerate states is reduced and the tunneling frequency of the 
system increases for the duration of the pulse. 

Furthermore, it is possible to tune the tunneling frequency of the qubit by 
applying a steady state current in the same manner as that of the a x operation. This 
allows a tuning of the quantum overlap of the two degenerate ground states in the 
qubit to a desired range. This is useful in an array of qubits where the tunneling 
frequencies of some or all of the qubits vary. By tuning the tunneling frequency of 
each qubit, the array can be tuned into a uniform range of tunneling frequencies, 
thus allowing more predictable application of quantum algorithms in the array. 

Tuning can be achieved for example on qubit 100 of Figure 8 by 
application of a steady state current to terminal 1 10-2 and allowing the current to 
escape through terminals 110-1 and 110-3. A method for tuning the tunneling 
frequency of a qubit would comprise controller 800 applying a steady state current 
at terminal 110-2 of junction 120, while grounding or allowing escape from the 
adjacent terminals 110-1 and 110-3. By varying the magnitude of the applied 
steady state current, the tunneling frequency of the qubit can be manipulated. 

Figure 9 shows a qubit array 900 including qubits 100-1 and 100-2. Array 
900 can include any number of qubits. Qubits 100-1 and 100-2, for exemplary 
purposes only, are each five-terminal qubits 100 as shown in Figure 8. Qubits 100- 
1 and 100-2 are entangled by junction 145 coupled between superconducting loop 
122-1 of qubit 100-1 and superconducting loop 122-2 of qubit 100-2. In order that 


qubits 100-1 and 100-2 can behave independently, it is desirable that junction 145 
can be opened to isolate qubit 100-1 from qubit 100-2 (and hence from all the other 
qubits in qubit array 900) and closed to entangle qubit 100-1 with 100-2, see 
Appendix A. Junction 145 can be controlled by controller 800 by applying a 
voltage to junction 145, for example capacitively through plate 146. 

Figures 10a and 10b show operation of plate 146 in opening and closing 
junction 145 and thereby entangling or isolating qubits 100-1 and 100-2 of array 
900. Controller 800 is electrically coupled to plate 146 at a terminal so that 
controller 800 can apply a voltage to plate 146, which is capacitively coupled to 
junction 145. Junction 145 can be any superconducting junction, including a two- 
dimensional electron gas, tunneling junction, or constriction junctions. 

As shown in Figure 10a, when controller 800 does not apply a voltage to 
plate 146, current can freely flow between superconducting loop 122-1 of qubit 
100-1 and superconducting loop 122-2 of qubit 100-2. Therefore, qubits 100-1 and 
100-2 are entangled. 

As shown in Figure 10b, when controller 800 applies a voltage to plate 
146, electrons are prevented from flowing through junction 145 by electric fields 
143 and junction 145 effectively becomes opened, thus isolating qubit 100-1 from 
qubit 100-2 . 

Figure 1 1 illustrates a qubit array 900 of five terminal qubits, where each of 
qubits 100-1 through 100-N is coupled by ajunction 145-1 through 145-(N-1) to 
its nearest neighbors. That is, qubit 100-1 through 100-N can be entangled by 
junction 145-1 and the entanglement can be switched by plate 146-1. Further, 
qubit 100-(N-1) and qubit 100-N can be entangled by junction 145-(N-1) and the 
entanglement can be switched by plate 146-(N-1 

Although array 900 of Figures 9 and 1 1 illustrate entanglements between 
qubits in a linear array, one skilled in the art should recognize that a two- 
dimensional array of qubits 100 can be entangled in this fashion. Each 
superconducting loop 122 can be switchably entangled to any number of other 
superconducting loops through ajunction 145 with a plate 146. 


Figure 12 shows a particular embodiment of qubit array with qubits 100-1 
through 100-N coupled by junctions 145-1 through 145-(N-1) and switched with 
plates 146-1 through 146-(N-1), respectively. In Figure 12, junctions 145-1 
through 145-(N-1) are each two-dimensional electron gas junctions. Additionally, 
junction 145-1, for example, is coupled such that it is within both superconducting 
loop 121-1 and 121-2, rather than simply linking superconducting loops 121-1 
through 121-2 as is shown in Figure 1 1. In this fashion, junction 145-1, through 
plate 146-1, controls the coupling of superconducting loop 121-1 with 
multiterminal junction 120-1, controls the coupling of superconducting loop 121-2 
with multiterminal junction 120-2, and controls the entanglement between qubits 
100-1 and 100-2. To avoid any interaction of the entangling junction with the flux 
in the qubit, the entangling junctions 145-1 through 145-(N-1) are isolated from the 
main loop of the qubit by a distance W145. This further decouples the qubit from 
any interaction with the environment or surrounding fields, thus decreasing the 
decoherence in the qubit. 

In operation, qubits and qubit arrays according to the present invention 
(such as, for example, the embodiments discussed above) are cooled to a 
temperature well below the superconducting transition temperature T c of all 
superconducting materials utilized to fabricate the particular structures. In an 
exemplary embodiment, the structures described in figures 1A through 13 are 
cooled to an operating temperature of about 10 milliKelvin so that all structures are 
superconducting, all phase shifters are operative and decoherence processes due to 
thermal fluctuations and inelastic scattering are suppressed. 

In accordance with current theoretical descriptions, for example, the 
Eliashberg theory of superconductivity (see, e.g., R. de Bruyn Ouboter and A. N. 
Omelyanchouk, "Macroscopic Quantum Interference Effects in Superconducting 
Multiterminal Structures," Superlattices and Microstructures, Vol. 25, No. 5/6 
(1999)) an order parameter describes current flow in superconductors and phase 
differences in multi-terminal junctions. Multi-terminal qubit 100 with intrinsic 
phase shifter 123 as shown, for example, in Figures 3 A through 3C have 
degenerate ground states of the qubit quantum system, the supercurrent circulating 


in superconducting loop 122 of qubit 100 being twice degenerate if no external 
magnetic field or transport current from any external current source is applied. 
The two degenerate states having the ground state energy and definite magnetic 
moment, |0> and |1>, correspond to minimal supercurrents circulating through 
superconducting loop 122 in clockwise and counter-clockwise senses. The two 
states associated with the supercurrent in superconducting loop 122 permit 
quantum computing in the standard fashion, which is described in many papers and 
books {see, e.g., J. E. Mooij, T. P. Orlando, L. Levitov, Lin Tian, Caspar H. van 
der Wal, and Seth Lloyd, "Josephson Persistent-Current Qubit," Science 285, 
1036-39(1999)). 

The role of phase shifter 123 in multi-terminal qubit 100 is to cause the two 
basis states of the qubit to be naturally degenerate. This is a major advantage over 
other qubit designs, for example that of Mooij et al., Science 285 1036, where it is 
necessary to apply a magnetic field in order to bring the basis states of the qubit 
quantum system into resonance (i.e., cause them to be degenerate). The magnetic 
field required by the system of Mooij et al., Science 285 1036, needs to be 
extremely finely tuned in order to maintain the resonance condition. This is 
because the precision to which the external field has to be tuned is approximately 
the tunneling amplitude A T between qubit basis states which is usually about 5 
GHz and corresponds to a magnetic field precision of one part in about 10 6 . 

One application of embodiments of multi-terminal qubits according to the 
present invention is a quantum computational random number generator. As a 
random number generator, the quantum states of an array of qubit 610 as shown in 
Figure 6A or array 680 as shown in Figure 6B or of array 900 as shown in Figures 
9 and 1 1 evolve to a state where the qubit quantum system on each of individual 
qubits (e.g, qubits 600-1 through 600-(N+l) of array 610) has an equal (or at least 
known) probability of evolving to each of the basis states |0> and |1>. The basis 
states, which are related to the superconducting current directions in the 
superconducting loops, are then determined, for example, by observing each of 
qubits with a magnetic force microscope, a SQUID magnetometer, a scanning Hall 
probe, or other magnetic probes, or alternately by applying transport currents and 


measuring voltage drops with, for example, apparatus 700 as described in Fig. 7, 
which will fluctuate if the state to be measured is in the higher energy state or 
remain static otherwise. Each determined state (clockwise current or 
counterclockwise current) corresponds to a bit value (0 or 1) so that the collection 
5 of determined states provides a random binary value having as many bits as there 
are qubits in the array. Quantum theory indicates that, with known (including 
zero) entanglements between individual qubits, a series of bits thus generated can 
be random without correlation or repetition between bits. 

Qubits according to embodiments of aspects of the current invention may 

1 0 also alternatively be read by other readout devices such as a magnetic force 
microscope (MFM) tip, a superconducting quantum interferometer device 
(SQUID) loop, or a Hall probe device. The readout device measures the weak 
local magnetic fields that the spontaneous supercurrents (clockwise or 
counterclockwise) cause in the vicinity of the multi-terminal qubit. More 

15 particularly, the MFM scans a microscopic magnetized tip attached to a cantilever 
across the surface and measures deformation of the cantilever as the mechanical 
force that acts on the magnetized tip. Alternatively, a superconducting loop can 
detect the magnetic flux in the vicinity of the multi-terminal qubit. Alternatively, a 
Hall probe can detect the magnetic flux in the vicinity of the multi-terminal qubit. 

2 0 Another possible read out system may use a difference in the absorption of 

circularly polarized microwave radiation due to the clockwise or counterclockwise 
currents by a multi-terminal qubit. 

The time required for a calculation and the interpretation of the read out 
results depends on the calculation performed. Such issues are the subject of many 

2 5 papers on quantum computing, for example P. Shor, "Introduction to Quantum 

Algorithms," Los Alamos preprint server condmat/005003 (April 29, 2000). The 
structures described herein can serve as components of quantum computing 
systems and also can implement any single qubit algorithm. 

In general a controller can be electrically coupled to provide current to the 

3 0 terminals of each qubit 100 in an array of qubits (for example, controller 800 

shown in Figures 8, 9 and 11). Controller 800 controls the currents in each of the 
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terminals and thereby is capable of controlling the initial states of each qubit, the 
entanglements between qubits, the application of magnetic fields to each qubit, and 
the measurements of the state of the qubit. As such, one skilled in the art will 
recognize that controller 800 can be a microprocessor based system operating 
5 software which controls the programming and readout of qubits single or in arrays 
as shown in Figures 6A, 6B, 6C, 9, and 1 1 . 

Although the invention has been described with reference to particular 
embodiments, the description is exemplary only and should not be considered 
limiting. One skilled in the art may recognize several obvious variations, which 
10 are intended to be within the scope and spirit of the present disclosure. One skilled 
in the art will recognize embodiments of other qubits according to the present 
invention which are within the scope of this disclosure. Various adaptations and 
combinations of features of the embodiments disclosed are within the scope of the 
invention. As such, the invention is limited only by the following claims. 
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We propose a novel superconducting flux qubit design based on mesoscopic 4-termmal junctions. 
The logical states correspond to different orientations of persistent currents in superconducting loops. 
Magnitude of the corresponding fluxes is of the order of a fraction of the flux quantum minimizing 
considerably interactions with the environment. No external fluxes are necessary to attain bistability 
or to perform logical operations. 


I. INTRODUCTION 

In the past twenty years, clear evidences have been gained that quantum information processing could offers signifi- 
cant advantages over classical information processing [1,2]. In parallel, it was recognized that superconducting systems, 
and particularly SQUIDs, are good candidates for the observations of quantum phenomenas at the macroscopic level 
[3]. A great body of experimental evidences has been accumulated to support this assertion. 

In the growing efforts to fabricate a quantum information processing device, SQUIDs are then natural candidates. 
Many proposals using either superconductor's phase or charge degree of freedom has been published, see [4] for a 
review. 

In this paper, a novel qubit design using the phase degree of freedom is introduce. Similarly to the rf-SQUID [4] or 
to the design introduced in [5] the logical states are represented by persistent currents of different orientations flowing 
in a superconducting loop. In the present design logical operations are not performed by manipulations of external 
fluxes. This is, among others, an advantage of this design. 

The paper is organized as follows. In Sec. II we present the layout of the multi-terminal phase qubit. In Sec. Ill 
the quantum dynamics of the flux degree of freedom in this device is studied and it is shown how it can be reduced 
to a two- level system (TLS). In the following sections, we go on showing how this TLS satisfies the requirements for 
quantum computation [6] : initialization and measurement in Sec. IV, a universal set of gates is given in Sec. V, VI 
and VII and finally estimates for decoherence are presented in Sec. VIII. 


II. MESOSCOPIC 4- TERMINAL QUBIT 

The basic constituent of the novel qubit design studied in this paper is the 4-terminal Josephson junction (Fig. 1) 
which has been the focus of extensive studies [7]. Nonlocal coherent transport in this structure was investigated in 
the companion paper (thereon referred to as Article I). The mesocopic 4-terminal junctions of figure 1 b) and c) on 
which we will focus our attention in this paper consist of four s-wave superconducting banks weakly coupled by a 
ballistic two-dimensional electron gas (2DEG) layer. As seen in Article I, in such devices the total current, I,-, flowing 
into the i-th terminal depends on the superconducting phases tp 7 at all the banks and has the form [8] 



where A 0 is the superconducting gap in the banks, the 7^ are geometry dependent coupling coefficients between the 
i th and j th terminals and fa = 0. Pair of terminals are either current biased or closed by a superconducting loop (see 
below). As noted on figure 1, we label the terminals such that terminals 1 and 2 sustain current I while terminal 3 
and 4 current J; below I and J will respectively be transport and persistent currents. 

As discussed in Article I, coupling between the I and J circuits is characterized by the matrix 

1:) » 
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(a) 




FIG. 1. 4-termmal Josephson junctions, (a) Junction with microbridges. (b) Mesoscopic junction with two-dimensional 
electron gas (2DEG) and symmetric current configuration. In this configuration, coupling between currents / and J inside 
the normal region is local (S = 0). (c) Mesoscopic junction in the asymmetric currents configuration. Coupling is non-local 


For multi-terminal junctions such that det( 7coup ) = 0, coupling between the currents J and J is local. Constriction 
junctions (Fig. la) and mesoscopic junctions in the crossed currents configuration (Fig. 2b) belongs to this category. 
In the later case : although coupling is really nonlocal, it is the extra symmetry in the current configuration that 
imposes det( 7coup ) = 0. Alternatively, when det( Tcoup ) ^ 0 coupling is nonlocal. This is realized for mesoscopic 
junction m the parallel current configuration (Fig. lc). In this case, and as emphasized in Article I, the phase across 
the two terminals in one current loop is affected by the direction of the current in the other loop. This will be referred 
to as phase dragging effect. As will be shown below, this effect plays a crucial role in manipulation, initialization and 
readout of qubit's states. 

B> connecting terminals 3 and 4 via a superconducting ring, a 4-terminal SQUID is realized. Figure 2 (the tt- 
junction will be introduced below) shows such a structure where, to take advantage of phase dragging effect, the 
asymmetric currents configuration of Fig. lc is used. 

In this configuration, the 4-terminal SQUID is reminiscent of the much studied rf-SQUID [9]. It is know that for 
rings of large enough inductance and under $ 0 /2 externally applied flux ($ 0 = h/2e, the flux quantum) the rf-SQUID 
can be bistable [4]. The bistable states correspond to different orientations of persistent currents or equivalently to 
different magnetic flux enclosed by the ring. Quantum behavior of the flux degree of freedom in this system as been 
the focus of numerous experimental and theoretical investigations over the past twenty years, see [4] for reviews. 

Bistability of the 4-terminal SQUID as also been the focus of theoretical investigations and as been confirmed 
experimentally [10]. For quantum computing purposes, a major advantage of the 4-terminal SQUID over the standard 
rf-SQUID is that bistability can be attained for arbitrarily small loop inductance. As a result, the enclosed flux is 
much smaller than for the rf-SQUID and weaker coupling to the environment is to be expected. Small inductances 
and thus small fluxes where the main advantage of the 3 junctions design of reference [5]. 

As in the rf-SQUID or 3 junctions case, an external flux of $ e = $ 0 /2 is needed to get bistability (see appendix 
A.). However, as will become clear in the next sections, this control parameter is not exploited to manipulate the flux 
(i.e. qubit) state. The external flux can thus be fixed to # 0 /2 for all qubits. This opens the possibility of replacing, 
as shown in figure 2, the external fluxes by a 7r-phase shifter in each qubit's superconducting ring. The net effect is' 
the same but this has the advantage that the 7r-phase shifter does not bring in extra coupling to the electromagnetic 
environment. 

Let us now proceed with the description of the 4-terminal SQUID. For clarity of presentation, the bulk of the 
calculations is given in the appendix and only the main results are presented here. It is useful to introduce the new 
variables: 


8 = ?*-<Pu <P = <P3~f4, X = 7;(<Pl+V2)-^{<Pi + 'p3), (3) 

where <j> is related to the flux $ threading the ring through <j> = 2tt$/$o mod27r. It is also useful to introduce the 
variable ip = $ + tt which, physically, is the phase related to the self-generated flux on the ring in the presence of a 
7r- phase shifter or, equivalently, a $ 0 /2 external flux. 
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k- Junction 


/ 


FIG. 2. Mesoscopic 4-terminal SQUID qubit. The 7r-phase shifter in the flux loop is added to attain bistability without 
external flux. The externally controlled transport current I affects the current J in the superconducting loop through the 
phase dragging effect and in turns the flux in the loop. 

In mesoscopic four terminal junctions where the 2DEG layer as a square geometry, the coefficient 7y are such that 
7i3 = 724 and 714 = 723. As shown in the appendix, this constrains x to take only one of two possible values, 0 and 
tc. We will thereon work under this condition but results can be generalized to the more general case. We define the 
new variables 7 = (7 13 + 7i4)/7i2 and 5 = (713 — 7i4)/7i2- The difference between local and nonlocal couplings lies 
in the value of 5, namely 5 - 0 for the local case but 6 ^ 0 in the nonlocal case. 

In appendix A, the free energy of the 4-terminal SQUID is derived for T = 0 and near the superconducting critical 
temperature T c . It turns out that the free energy has the same form in both cases with differences only in the values 
of their coefficients. We therefore expect the same behavior to hold at all intermediate temperatures as well. All 
discussions in the main body of this article are devoted to the T = 0 limit keeping in mind that the arguments are 
more general and that the same expressions can be derived near T c as well. Evidently, the limit T — 0 is more 
pertinent for qubit operation where one should expect longer coherence times. 

i,From (A18) (keeping only the cp dependent terms), the free energy for the flux loop can be written as 


where C— 2irI 0 L/$o is the dimensionless inductance with I 0 the Josephson critical current between terminals 1 and 
2 (see appendix A). Energy is expressed in units of the Josephson energy E 0 = hl 0 /2e of the 1-2 loop and current in 
units of io- The coefficient of the second term, £(J"), and <po(I) are functions of the transport current I : 


with ¥>o(0) = 7£/v / 2 and e 0 = 26. The energy (4) has two local minima at tp = ±<po{I) with energy difference 2e(J). 
These minima correspond to different orientations of the persistent current in the superconducting loop. It is clear 
from (4) that the system can be bistable even if the superconducting loop as a very small inductance L. As stressed 
earlier, small inductance results in a small separation between the minima and correspondingly a smaller coupling to 
the environment. 

The main difference between locally and non-locally coupled 4-terminal junctions, for which 5 = 0 and ^ 0 respec- 
tively, now becomes apparent. In the former case, one as eo = 0 and the minima are always degenerate. In this case 
(or when e{I) is tuned to zero) the barrier height between the minima is given by 


In the non-locally coupled case, the bias term can be used to breaks this degeneracy and make one of the minima more 
favorable than the other. Hence, bias energy and barrier height can be tuned by the externally applied transport 
current. In the following sections, we will show how these properties can be used for initialization, readout and 
manipulation of the qubit's state. 


UM) = ^ (M - Vo(/)) 2 - e{I) sign(^), 


(4) 




(6) 
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III. QUANTUM DYNAMICS 

In the presence of a small capacitance C across the junction 1-3, the system's total energy is 

E=^ + U,(I), (7) 

where Q is the number of charges on the capacitor and C = E 0 Cf{2e) 2 the dimensionless capacitance. The system's 
Hamiltonian is obtained by replacing Q by the appropriate operator in (^-representation, Q -> id /dip. In units of E 0 
one then obtains 

H = S + h (M ~ Vo(/))2 " e < / ) si S n (v)- (8) 

This is equivalent to a particle of mass proportional to the capacitance C moving in the double well potential of 
equation (4) . For small enough capacitance, this system can behave quantum mechanically. 
Near the two minima of (4), the classical oscillation frequency, in units of Eo/h, is given by 


' %/£C EqVLC' 


0) 


If this frequency is much smaller than the barrier height, in other words if w 0 <C £/&(/), then the low energy levels on 
both sides of the barrier can be approximated by 


E n , ± = + w 0 (n + 0 . (10) 

leling between the r 
I by the tunneling sp 

( "oCvo{I) V'\_ A 


When the minima are degenerate, coherent tunneling between the right and left well is allowed. In this situation, 
degeneracy between the two lowest levels is lifted by the tunneling splitting. 


So = dtpy/2CU v (I), A — p + c 
For the potential (4), A is identically zero and 


? I" coqC 1 j 


One then obtains the tunneling splitting 

A(J) = 2w 0 QcW)V2c) ^ 2 e- 2U >W'»°. (14) 

In the regime were the attempt frequency uiq is much smaller than the barrier height, but much larger than tunneling 
splitting, energy bias e(T) and temperature, truncation to the Hilbert space spanned by the system's two lowest energy 
levels is possible. In the basis of the ground states of the right and left well respectively (which will be denoted {|0), jl)} 
for obvious reasons), this yields the usual two level system Hamiltonian 

a,* = A{I)a x - e{I)o z . (15) 

The u's are Pauli matrices and the coefficients and A(I) are given by (5) and (14) respectively. We have thus 
mapped the 4-terminal SQUID flux state's to a pseudospin. We will now proceed to show how this pseudospin can 
be used as a qubit. 
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IV. INITIALIZATION AND READOUT 

Initialization is the process of taking an arbitrary state and bringing it to a well-defined state (which is usually 
chosen to be |0)). This evidently is a non-unitary process and is therefore possible only through coupling to an 
environment. Obviously, measurement is also a non-unitary process and also requires coupling to an environment (the 
measuring apparatus). 

In this design, initialization and readout are performed through the phase dragging effect by applying a transport 
current. More specifically, for initialization a transport current / is applied for sufficiently long time. Because of the 
resulting energy bias (5), the system will relax to the lower energy state after sufficiently long time. By choosing the 
direction of current one can choose to whichever state to initialize to. 

To readout the state of a qubit, one can use the fact that the critical current in the transport current loop depends 
on the state of the flux loop. As derived in appendix A, for the system to be bistable, 9 (denned in (3)) must satisfy 
the constrain \9\ < & max , where ff max is given by (All). When this is satisfied, the current-phase relation in the 
transport current loop is given by (A15) and is reproduced here for convenience 

T c- 6 5 0 ■ , x 

J = Sln -__ cos _ slgnM (16) 

The critical current therefore occurs near 6 = tt with the value I c = 1 T 6/2, depending on the qubit's state. In 
order to measure this state, one applies a current 7=1. If the qubit is in the state + Vo , the current exceeds the 
critical current and the system goes to the non-stationary state with a voltage drop across the transport terminals. 
Un the other hand if the system is the state -<p 0 , the current will be below the critical current and no voltage will 
be detected. b 

An important condition here is that when & « w the system should still preserve its bistable behavior. Therefore 
we need to have 6 max > tt. This imposes a minimum value for the equilibrium flux <p o (0) > 4^2(5 H), or equivalent^ 
puts an upper limit on S - 

5 < £ 7 2 /8. (17 ) 

Minimization of the flux ^ 0 while respecting (17) requires 8 < 7 or equivalently 7l3 a lu . 

Note that initialization can alternatively be performed by measuring the qubit's state. If the measurement result 
is jO), initialization is complete. Otherwise, one can apply a it pulse (see below) to obtain jO) (this can be followed by 
a second measurement to cope with any imperfections in the pulse application). 


V. SINGLE QUBIT OPERATIONS 

Single qubit operations are generated by the single qubit Hamiltonian (15). As stressed in sections II and III the 
coefficients s{I) and A(7) of this effective Hamiltonian can be controlled through the transport current 7. This gives 
the possibility to coherently manipulate the qubit's state. 

Let us first consider locally-coupled multi-terminal qubits for which 5 = 0. In this situation, e(7) is fixed to zero and 
only A(7) is externally controllable. For coherent manipulations, the control current J should not exceed the critical 
value 7 C a 1 where the system becomes non-stationary. In practice, the maximum allowed value for 7 is slightly less 
than this critical value in order to avoid tunneling and thermally activated processes in ^-direction. Hence, currents 
roughly in the range [0, 1[ are useful for coherent control of A (7). In this range, the barrier height U b (I) changes by 
a factor of 1/2. Since tunneling amplitude A(7) is exponentially dependent on this parameter, reducing U b (I) by a 
factor of 1/2 can have a significant effect on the value of A(7). Thus, one can practically switch tunneling on and off 
by applying transport currents 7 = 1 and I = 0 respectively. As a result, when 7 = 0, the qubit does not evolve, it 
is frozen . With a transport current 7 = 1 however, the effective Hamiltonian is A(1W S and the qubit's evolution 
operator is 

X{a) = e~^l\ (18) 

with a = 2E 0 A(l) t/h and where we reintroduced proper units for the tunneling splitting. Choosing t and 7 such 
that ol tt, one obtains the x pulse of the previous section. Therefore, for 6 = 0, a qubit is either frozen (7 = 0) or 
undergoes coherent oscillations (7 « 1). 

K . In ^ asymmetric current configuration of figure 2, <5 ^ 0 and the transport current will also produce an energy 
bias. When the bias exceeds the level width, tunneling between the two basis states stops because they are out of 
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degeneracy [11]. Hence, if I = 0, the levels are degenerate but the qubit is frozen because the barrier height is at it's 
maximum. When 7=1 however, e(I) reaches it's maximum value V26 and the barrier height it's minimum value. 
Coherent tunneling is nevertheless prohibited because the degeneracy condition is not satisfied. As a result, for 6 ^ 0 
a qubit can either be frozen (7 = 0) or the effective Hamiltonian E{I)cr x is generating a phase difference between the 
basis states (I > 0) 

Zifi) = e-*'*' 2 , (19) 

with $ = 2£ 0 e(7) t/h. 

It is well known that to implement arbitrary single qubit logical operations, one needs rotations around at least two 
orthogonal axis [12]. Since the situations 6 = 0 and ^ 0 cannot be realized in the same 4-terminal SQUID we clearly 
not yet have reached this goal. Different ways to remedy to this are possible. First, for 6 ^ 0 the 4-terminal SQUIDs 
can be engineered such that for 7 = 0 tunneling is not suppressed (e.g. by working with smaller capacitances). In 
this case, both X(a) and Z(/3) operations are possible. It was shown in [13] how, in that situation, to 'freeze' a qubit 
evolution using NMR-like refocusing sequences. Alternatively, for 8 = 0, it will become clear in section VII that by 
using two-qubit gates and the single-qubit gate X(a), one can implement arbitrary single-qubit gates. Finally, we will 
see in the next section how the 4-terminal SQUID can be modified to effectively combine both the 6 = 0 and <5 ^ 0 
cases. 

Before moving on, some comments are in order. First, for <5 = 0, when setting the transport current to 7 = 0, one 
would expect the qubit to be frozen. However, since the tunneling amplitude A(0) is not identically zero, care must 
be taken to prevent residual coherent tunneling from introducing bit flip errors. Accordingly, we require that the time 
for a bit flip when the qubit should be frozen hir/2E 0 &.(0) be much larger than the time 7i7r/2.EoA(l) when coherent 
oscillation are desired (when I k 1). For the error probability to satisfy the threshold condition for fault-tolerant 
quantum computation [2], one then roughly requires that A(0) < 1(T 4 A(1) or equivalently that U b (0) > 41nl0w 0 , 
which in turns puts a limit on the size of the capacitor 

„ (4 In 10) 2 

c> mw- (20) 

Moreover, it is important to make sure that the bias energy is always less than the difference between the energy of 
the lowest energy states and the first exited state. Otherwise, there might be leakage out of the computational Hilbert 
space causing loss of information [14]. To avoid this undesirable tunneling we require 5 < w 0 /\/2. 

Finally, the double well potential (4) is obtained from (A18) which is a function of ip and 6. It is through relation 
(16) that 6, and thus U^, is a function of the transport current 7. To obtain (4) however, the small term proportional to 
5 in (16) was dropped. For very small 5, this is reasonable but care should be taken because, through this term, 0 does 
not only depend on the transport current I but also on the qubit's state. Since 6 is not an observable quantity, this 
will not lead to additional environmental coupling. However, this effect, although very small can cause computational 
errors after a large number of operations. To minimize this, the junction parameter 6 should then be chosen as small 
as possible. Since we are interested in junctions such that 6 <C 7 <IC 1, this corresponds to 713,714 <U 712 and is not 
hard to achieve in practice. Note however that while this later requirement is in agreement with (17), junctions with 
very small <Ts will require a correspondingly high level of accuracy in the applied transport current for measurement. 


VI. MESOSCOPIC 5-TERMINAL QUBITS 

As exposed in the previous section, junctions with 6 = 0 and 5^0 offer complementary possibilities for single bit 
operations and it would be practical to combine those possibilities in a single qubit design (this would for example 
remove the need of refocusing to freeze qubits in the asymmetric current configuration). This is impossible for the 
single 4-terminal SQUID because, as seen in figure 1 b) and c) , non-local and local coupling of currents inside the 
2DEG layer are realized in differently engineered junctions. 

We suggest here to use the 5-terminal junctions illustrated in figure 3. Coupling between the two flux rings is used 
to create entanglement and will be discussed in the next section. The extra lead for transport current in this design 
gives the opportunity to use either symmetric (6 = 0) or asymmetric (6 0) current configurations. More specifically, 
the qubit on the left hand side of figure 3 (where the middle junction is disconnected from the current source) is 
completely equivalent to the asymmetric 4-terminal SQUID. In this configuration, the transport current tunes e(J) 
(as discussed in the previous section, it will also tune A (7) but as long as the truncation conditions hold, this does 
not change the dynamics). On the other hand, the rightmost qubit in figure 3 is in the symmetric configuration 
and therefore 6 = 0. For the current values indicated on the leads of this qubit, the description is analogue to the 
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FIG. 3. Two coupled five-terminal qubits. The left qubit has asymmetric current configuration while the right one has 
symmetric one. 



FIG. 4. System of N coupled qubits. The distance bet-ween the flux regions is maximized to have less magnetic interaction 
between the qubits. 


4-terminal case already discussed. Therefore, in this symmetric configuration, the transport current tunes A (7) while 
e{I) is fixed to zero. 

As a result, without transport current, the 5-terminal qubit is frozen (with appropriately chosen junction parame- 
ters). Applying a transport current I as 1 in the symmetric configuration a A' (a) logical operation is performed while 
the asymmetric configuration yields a Z{p) operation. Arbitrary single qubit gates are then readily implementable in 
this design. This is done without manipulation of external fluxes but by injecting currents in appropriate leads. 

Note finally that measurement and initialization should be performed here using the asymmetric current configu- 
ration. 


VII. TWO-QUBITS OPERATIONS 

To obtain a universal set of gates for the 5-terminal qubits, a prescription for a non-trivial qubit-qubit interaction 
is needed [12]. Here again, this interaction is provided by the phase dragging effect. Two qubits are coupled through 
an additional mesoscopic 4-terminal junctions as presented schematically in Fig. 3. 

The interaction can be turned off by applying a gate voltage to the 2DEG region between the two qubits. This 
voltage depletes the middle region and, as a result, prevents coherent coupling between the two flux loops. This 
interaction is studied in more details in appendix B and the resulting effective two-qubit Hamiltonian is 

H 2 gb = H!+H 2 +n <7l,z ® d2,2 , (21) 
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where i = 1,2 indexes the qubits, H, is the one-qubit Hamiltonian (15) and 



(22) 


The primed coefficients axe parameters of the connecting 4- terminal junction (see appendix B). As indicated in the 
previous sections, the single qubits Hamiltonians can be manipulated individually through their transport current 
leads. When no transport current is applied on either qubits and when the 2DEG connecting the two qubits is not 
depleted, the system's evolution operator is then given by 


where £ = 2E 0 £lt/h and CP stands for 'Conditional Phase shift'. 

As schematically shown in figure 4, N qubits can be connected in this way to form an iV-qubit quantum register. 
From the 2-qubit logical operation CP(Q and the single-qubit gates defined earlier, it is possible to implement arbitrary 
A r -qubit logical operations on such a register [13]. For example, a Controlled-NOT gate can be implemented, up to 
an irrelevant global phase factor, by the following sequence of one- and two-qubit gates 


Some comments are in order. First, magnetic interaction between qubits should be very small to eliminate de- 
coherence due to unwanted qubit-qubit interactions. The qubits should therefore be disposed in such a way as to 
keep the regions containing the maximum flux farthest apart. Secondly, it is possible, in this design, to apply logical 
operations simultaneously on distinct qubits. Since it is known that only quantum computers having this kind of 
classical parallelism can hope to satisfy a threshold condition for fault-tolerant quantum computation [2,15], this is a 
major advantage of this design. Note that not all (solid-state) quantum computer architectures have this ability [4]. 

Moreover, figure 4 shows an A'-qubit array in it simplest design : a linear array of qubits limited to nearest-neighbor 
couplings. While more complex coupling schemes are possible, this configuration is, at least for some computational 
tasks, sufficient [16]. Finally, as announced in section V, for two-coupled 4-terminal SQUIDs in the symmetric 
configuration (i.e. where 6 = 0 for both qubits), the two-qubit gate CP(Q and the single-qubit gate X(a) are 
universal. This is quite obvious because CP(Q gives the possibility of accumulating relative phases, something which 
was absent when only X (a) was available. 


As for any other qubit design, the most important limitation for performance is decoherence due to interaction of 
the qubits with the environment. Many sources of decoherence that may be active for superconducting flux qubits 
were discussed in [17], see also [4]. These are fluctuation of background charges, quasiparticule tunneling in the 
superconductors, coupling of fluxes to nuclear spins, electromagnetic radiation and finally unwanted magnetic qubit- 
qubit couplings. At low temperature and for small equilibrium fluxes, coherence time was shown to be limited by the 
later source [17]. In the present design, this can be minimized as was discussed in the previous section. 

Let us now consider additional sources that may cause decoherence in this system. One potential source of de- 
coherence is the fluctuations of the external transport currents. These cause fluctuations of the tunneling matrix 
element A(I) and of the energy bias and in turn affects the relaxation and dephasing time of the qubit [4,20]. 
The effects of fluctuations of a general external control field on the relaxation and dephasing times of a qubit are 
described in [4]. We will consider here fluctuations of transport current in the asymmetric configuration where they 
should affect dephasing and relaxation more strongly. In this situation, fluctuations of currents 51 will provide a term 
of the form <r z 5I in the system + environment Hamiltonian. The power spectrum of the current fluctuations is given 
by Johnson-Nyquist relation 


(23) 


CN 12 = X 2 {Zir/2) CF 12 (tt/2) Z 2 (tt/2) X 2 (k/2) CP 12 (ir/2) Z 2 (it/2) Zi(ir/2). 
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VIII. DECOHERENCE 



(25) 


where Z(cj) is the impedance of the transport current circuit. For a purely real impedance, Rj, the strength of 
qubit-environment coupling is characterized by the dimensionless parameter [4,21] 
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where R K = h/e is the quantum resistance. Hence, to minimize coupling, J2, should be as large as possible and S 
small. This last constraint goes along the constraints already imposed on 6 in sections IV and V. From (26). one can 
estimate the relaxation and dephasing times respectively given by 


) coth ( ^ 


(28) 


where AE - y/e{I) 2 + A(I) 2 and T cs is the temperature at which the current source circuit operates, which is not 
necessarily equal to the working temperature of the qubits (T). T cs is usually higher than T. In the large T cs limit 
the two terms in (29) have the same functional form and we find 

The relevant parameter for quantum computation is less the dephasing time by itself than the quality factor Q ■ 
the ratio of the dephasing time over the characteristic time it takes to implement a single logical operation. For the 
asymmetric current configuration this characteristic time is h/E Q 8 and we obtain for the quality factor 


W nS{R K J {k B Tj ■ 


(30) 

For typical junction parameters, the Josephson energy E 0 = hl 0 /2e is of the order of XX. For the corresponding 
frequency of operation of the qubit, the lead impedance is of the order of the vacuum impedance and thus Rj ~ 100 fi 
d /» Smce f 3 ^ for workin S temperatures in the miliKelvin range, the quality factor should be large. The ratio 
R, JJr K ~ 10 1S the hmitlI1 S fact °r here. However, as we shall shortly see, this ratio could potentially be modulated. 

The finite conductance of the normal quasiparticles through the 2DEG can also cause decoherence. The coupling 
of the quasiparticle current to the system is through the term I^Mfyz- Following exactly the same procedure as 
above we find 

R K 

a= (2^7° (J )> (3D 

where i? qp is the quasiparticle resistance of the junction (or shunt resistance in the Resistively Shunted Junction 
(RbJ) model). At small temperatures {k B T <C AE) we find 

r-i-ir- 1 -VMllf 11 *} f AE \^{ AE \ <pl(I) { R K \ { AE\ 

T * ~ 2 r «- ~ 32^ [r-J {—) C ° th {j*T J K -id{RZ){—)- W 

Although (32) does not show explicit T dependence at low temperatures, R qp depends on temperature. As shown in 
[18] the quasiparticle resistance is Rq P ~ T cosh 2 (E s /2k B T), where ±E S are the energies of the Andreev bound states 
inside the normal region. R^p exponentially diverges as T 0. Thus by working at low temperatures, this source 
of decoherence is made negligible. This is along the same line of reasoning as the discussion in [17] on decoherence 
due to quasi-particules tunneling in tunnel junctions although E s can be much smaller than the gap A 0 , which is the 
energy scale that comes into the exponent in tunnel junctions. 

Moreover, care must be taken that the ac Josephson voltage generated by the phase fluctuations does not generate 
mter-Andreev level transitions [19]. This imposes u 0 < E s which therefore requires that C > 1/LE 2 . This constrain 
on the capacitance does not contradict the condition (20) already imposed in section V to minimize* bit flip errors. 

Pairs of qubits are also coupled to the electromagnetic environment through the gate voltage used to tune qubit- 
qubit interactions. When this interaction is turned off, i.e. when the coupling 4-terminal junction of figure 3 is 
depleted, small fluctuations of the voltage should have no effect on qubits as long as the gate voltage is large enough 
This is because the fluctuations of voltage introduce fluctuations in phase y?' (through the parameter see appendix 
-J 1106 Jf 5 V (K ' > theSe fiuctuation cou P le ver y weakly to the state of qubits and therefore do not introduce 
considerable decoherence. Moreover, these fluctuations can be minimized by working with a low resistance voltage 
source [4,20]. When the voltage is turned off, the region is no longer depleted and the qubit-qubit interaction is on. 
In this case, since the gate voltage is zero, the voltage source can be practically disconnected from the circuit and 
therefore no decoherence will be introduced. 
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An additional source of dephasing is the spread in qubit parameters and the errors in gate applications that follows. 
This is common to all solid state architectures and can be taken into account, for example, by calibration [24]. 

Decoherence in this design is then limited by fluctuations of transport currents. With the estimated working 
frequencies, about XX logical operations can be performed within the dephasing time due to this source. 

An advantage of this qubit design is that logical operations are not performed by manipulations of externally applied 
magnetic fluxes. For rf-SQUIDs and for the qubits introduced in [5], independently controlled fluxes of the order of 
$o/2 must be applied to two different regions of each qubit. Given the size of of those structures this is a technically 
challenging issue and imperfect shielding will cause dephasing. 

Moreover, in the design studied here, since the flux is not used to manipulate qubit's state it can be replaced by 
7r-phase shifters as discussed in section II. Inclusion of these elements can complicate the fabrication, as this will, for 
example, necessitate the use of d-wave superconductors. Alternatively, static external magnetic fluxes can be used. In 
this case however, fluctuations of the fluxes provide yet another coupling to the environment. As long as the spread 
in loop sizes is not to large, a single magnetic field can be applied to all qubits. In this case, the qubits are coupled 
to the same electromagnetic environment and, encoding two physical qubits in one logical qubit, the techniques of 
decoherence-free subspaces can be used [25,26]. These techniques where implemented successfully to protect quantum 
information from decoherence in situations where the qubit-environment coupling symmetry was the same as the one 
here [27,28]. 

Finally, one of the great challenges in the implementation of an useful quantum computing device is that coupling 
to the environment should be as weak as possible during coherent manipulations (i.e. computations) but should be as 
strong as possible during initialization. This later requirement comes, in part, from the fact that many error correction 
techniques requires continuous supply of fresh ancilla qubits. Therefore if initialization is slow (i.e. coupling to the 
environment is weak), the logical qubits for which error correction is intended will have had the time to decohere before 
fresh ancilla are prepared [6]. A possible solution is to initialize all the necessary qubits beforehand and distribute 
them in the quantum register as needed during computation [6]. This will however require very large quantum registers 
and many manipulations of qubit's state. Another possibility, is to tune the individual qubit's environment as needed 
either for computation or initialization. A possible way to implement this tunable environment is to use as transport 
current leads one-dimensional arrays of dc-SQUIDS as was done recently in [23]. By applying an external flux, the 
array act as a tunable electromagnetic environment where R T can be changed by many orders of magnitudes. Hence, 
during coherent manipulations using the asymmetric current configuration, one chooses a high resistance environment 
while during initialization, a small Rj. If this can be done on individual qubit, one can then tune the decoherence 
time of individual qubits by many orders of magnitude. 


We described in this paper novel superconducting flux qubits. These qubits are built from mesoscopic 4-terminal 
junctions. The logical states are represented by different orientations of persistent current in a superconduting loop. 
The corresponding fluxes differ only by a fraction of the flux quantum and therefore only weak coupling to the 
environment is to be expected. Logical operations are performed, through the phase dragging effect, by applications 
of transport currents. This compares favorably to other superconducting flux qubits [4,5,24], where two independent 
controllable external fluxes must be coupled to each qubits. 

Measurement is performed by applying a precisely tuned transport current in the control leads. As opposed to the 
other flux qubits, there is then no need to couple an additional dc-SQUID to the qubits and consequently simplifies 
fabrication and manipulations. 

Decoherence should be slow in this system allowing for a large number of logical operations to be applied. To 
help in further reducing the impact of qubit-environment couplings, encoding physical qubits into decoherence free 
subspaces [25,26] or the use of active error correction techniques [2,29] is possible. 


We start from the current-phase relationship (1) to derive the free energy for the four terminal SQUID. The 
expression (1) is simplified considerably near T = 0 and near T e . Near T = 0 one finds 


IX. SUMMARY AND CONCLUSION 


APPENDIX A: THE 4-TERMINAL SQUID 



(Al) 
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The Josephson coupling energy of the junction, Ej, which is related to the supercurrents U in Eq. (Al) through the 
relation J, = (2e/h)dEj/d<pi is given by 


Ej = 


Near T c on the other hand, one can write 


I t 


with the Josephson coupling energy 


Ej = 


In the subsections A and B we derive the free energy of a four terminal SQUID in these two limits. By comparing 
the results of the two subsections, one can easily see that the main properties'are the same in the two extreme cases, 
T « 0 and T « T c . Therefore one can conclude that similar behavior should be expected at all other temperatures. 
For simplicity of calculations, in this article more focus on temperatures close to zero, which is the appropriate limit 
for qubit operation. 

Considerable simplification in the calculation will be achieved if 

713 =724, 714 =723 (A5) 

In that case x, as defined in (3), can only take two values x = 0 or ir (see subsections A and B). This condition 
naturally happens in a mesoscopic 4-terminal junction with square geometry. 


-s^e^K 2 ^)! (A2) 


TTAo(r) 2 ^ . 


(A3) 


h tAq(T) 2 ^ 


1. The free energy near T = 0 

At T = 0 the Josephson energy of the 4-terminal junction will have the form 

= 2 ^ | cos ( l£L Y £i ) | (A6) 

The energy is in units of E 0 = hl 0 /2e, the Josephson energy of between terminals 1-2, and Jo = 7rAo/2e is the 
Josephson critical current between those terminals at T = 0. If we use the definitions of Eq. (3), the coupling part of 
the Josephson energy becomes 

*-> - - 2 M~ - I)h» h - f )l + - h + f)h>< h {*¥ + 1)1 ] 

We use the simplifying assumption (A5) to get 

£cou P = - 4 713 Max{|cos (*±*) cos || , |sin (*±£) sin f |> 
-4 7l4 Max{|c 0 s(^)cos§| , | sm (^sin§|} 
In the presence of a ir- junction or an external flux <£ e = it, we substitute <j> -+ <p = 4> — ir 

E m »p= -4 7ls Ma*{|cos(i±f±£)cosf| , |si„ (l±|±«) !in ||, 

— H-t^HI • K( =± f =£ )-»fl> 
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For small ip and 9 

«(— i— J-->(— f-J-7I <A7 ' 

In that case, the minimum of the energy with respect to x happens when x = 0 or jr. Otherwise x nas t0 De close to 
jt/2 and adds an additional undesirable factor of l/\/2. Therefore we can write 

= -iM^y* jcos (^±ft£) | + Tl . |cos (i±p») | , 

Tl 3j S i„(l±|±l)| +7l4 jsi„(l±p£)| } (A8 , 

Substituting j = j i3 + 7 14 and 5 = 713 — 714 and after some manipulations 

E C0lip = -2Max{( 7 + 6) |cos (^f^) | + (7 - f) |cos (l±p«) | , 

(7 + 6) jcos (iZjp?) j + (7 _ ,) | cos (,1Z|±£) | } ' (A9) 

When 6 = 0, the first (second) argument of Max is the largest for negative (positive) <p. Changing the sign of cp, the 
maximum switches from one argument to another. This results in two minima for the energy at two ip's with opposite 
signs. At finite 6, the two minima still exist as long as \6\ < 6 max , where 6 max is the value of Q at which the two 
terms inside Max become equal. 0 mai can be obtained by solving 

(A10) 


where <p(0 max ) is the equilibrium value of <p at $ = 6 max . Using the minimum value for equilibrium <p we find 

9 max = 4tan- (|tan^) « 4tan"> (^) (All) 
By making 5 <tC 7 we can make & max large enough so that we have bistability for all -tt < & < tt. We then can write 

= -4 7 cos j cos (?-=jM) - 4,5 sin £ sin (^jM) signal (A12) 
Assuming small tp, the free energy of a 4-terminal SQUID, can therefore be written as 

u = 27 " 19 ~ 2 co 4 " kM ~ 2V2j (* + ^) cos f " 2V2S (* ~ ¥) sin \ Sign( ^ (A13) 

We are interested in the limit where 6,k <C 7 <t; 1. Minimizing with respect to 6, we will have 

The second term adds a small correction to the critical current which can be neglected. The last term on the 
other hand, changes sign from one minimum to another. This is important in the readout of the final state. The 
current-phase relation therefore can be written approximately as 

\ sign(^) (A15) 
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The last term in (A15) can cause computation error when large number of operations is involved. Thus we want it to 
be small. Neglecting the last term we find 

^ = 2sin- 1 (7) ( A16) 
The ^dependent part of the free energy will then be 

Uv = 2C~ C ° S I " 2v/ ^ sin 5 si sn(v) (A17) 

We can write this energy as 

UM) = {M - VoCO] 2 - e{I) signfo>) - U b (I) (A18) 

with 

^(7) = ^!) ( 1+ yrrp) 1/2 ; £{ i) = s 0 i(i+vT^py' / \ i7»(j) = i£/»(o)(i+>/rr75) (A i 9 ) 

But now <^o(0) = -yC/V2, e 0 = 25 and £4(0) = <p 0 {0) 2 /2£. 

2. The free energy near T c 

Near T c , we simply write the Josephson energy as 

Ej = ~ T« cos ^' _ Pi) (A20) 

Here the energy is in units of E 0 = hl 0 /2e), the Josephson energy of between terminals 1-2, with 7 0 = tt7i2 Ag/4er e 
being the Josephson critical current between those terminals. We also normalize all other 7^- to 712- Introducing the 
new variables as in Eq.(3), we can write the Josephson energy of the junction in units of E 0 — ftj 0 /2e as 

Ej = - cose - 734 cos 4> + E coup (A21) 

where 

- - [„,«. (4±2 - I) - (iz£ _ I) +7UC05 (±z£ + 1) +toC0S (i+f + a) ] 

Assuming the condition (A5), and defining new parameters 7 = 713 + 714 and 8 = 713 - 7^ we find 

E coup = -2 [7 cos - cos - - J sin - sin - J cos x (A22) 

Minimization with respect to x leads to x = 0 or it depending on the sign of the expression inside the parenthesis. 
This is directly consequence of (A5). In cases more general than (A5), x can take any value (other than 0 and ir) and 
the expressions become more complicated. We therefore find for the total energy of the 4- terminal junction 

Ej = - cos0 - kcos</>- 2 ^cos B - cos ^ - <Ssin ~ sin | j (A23) 

with k = 734. The free energy of a SQUID made of this 4-terminal junction is given by 

U = ^ 2£^ ~ 19 ~ COS0 ~ KCOS< t > ~ 2 Jt cos \ cos I - * sin I sin (A24) 

Here, C = 2irI Q L/$ 0 = 2eI Q L/h is the dimensionless inductance of the ring. Now let's assume that an external 
magnetic flux <j> € = tt is applied to the SQUID loop. Changing the variable ->■ tp = <j> - it, will give 
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x - 2tan~ 


(A25) 

Exactly the same equation is obtained if one adds an ideal 7r-junction in the SQUID-loop instead of applying the 
external flux. In that case <p will correspond to the actual flux treading the loop {<p = 27r§ /$ 0 }, and 4> (= tp + w) will 
be the phase across the junction. 

The interesting regime for qubit considerations is n < S, 7 <§: 1 and -y£ < 1. In this limit, the effect of the transport 
current on the qubit is maximal and the effect of the qubit in the current is minimal. The last condition insures that 
the difference between the flux produced in the two states is small and the qubits are so called "quite" ; they do not 
interact with each other. The (^-dependent part of the free energy can then be written as 

V 2 „ I 0 10 8 i£>| 

u v = 2£ - 2 |7COS-sin-+5sin-cos| (A26) 

When the first term in the argument of the absolute value is the largest, the free energy remains unchanged when 
ip -> -ip. In that case, the free energy has two minima at two values of (p = ±<p 0 . Thus in order to achieve bistability 
the first term inside the absolute value should be the dominant one. This is possible only if |0| < 8 max where 

Notice the similarity between this equation and (All). At small <p we can write 

tp 2 I 8 8\ 
U * = 2£- UP cos 5 + 26 sin 2 ( A2S ) 

Assuming that the sign of the argument inside the absolute value is determined by the first term, Eq. (A28) can be 
approximated by 

U v = 2£~ 7M cos \ ~ 26 si S n (v) sin ^ (A29) 

We find two minima at y> ± = ±jCcos(8/2). The last term breaks the degeneracy between the two minima. The 
energies of the states become 

E± = - )p 2 L cos 2 9 - t 2<5 sin | (A30) 

The energy difference between the levels will be 2e{8), where 

£ (0)=2<5sin| (A31) 

This term can provide a z operation for quantum calculations. 

The phase 8 is determined by the transport current i". From Eq. (A25), the minimization dU/88 = 0 leads to 


The second term in (A14) is the same at both minima. The third term on the other hand, changes sign from one 
minimum to the other. This small correction is important in the readout of the qubit state but it causes error in 
quantum calculations. Since 7, 5 <C 1, we can approximately write 

0 = sra _1 (J) (A33) 
In terms of the current /, we can write the energy as 

U V V) = 2£ IM - ¥>o(/)f - e(i) signfoO - U b {I) (A34) 
ipo, £ and the barrier height Ub{I) depend on the transport current I through 

*<>( J ) = ^( 1 +v'l-J 3 ) 1/a , £(J) = £ 0 /(l+v / !^)~ 1/2 , U b {I) = \u b {0) (l + vTTf*) (A35) 

where ip o (0) = fC, s 0 = 25 and U b (0) = tp 0 (0) 2 /2C. These expressions are exactly the same form as what we found 
near T = 0. 


-69- 


APPENDIX A 


APPENDIX B: COUPLING TWO QUBITS 

Let's consider two qubits coupled via a mesoscopic 4-terminaI junction. Each qubit has its own control multi- 
terminal junction for manipulation and readout. The free energy of the coupled qubits is 

rT (fa + fa — 7r) 2 (fa + fay — 7r) 2 „ , 

U = 2£ + 2£ + Ej ^^ + E M<h,6 2 ) + E l2 (<f>[,fa 2 ) (Bl) 

Eji and E J2 are the Josephson energies of qubits 1 and 2 respectively. The qubits have the same parameters k, 7 
and 5. E n is the Josephson energy of the coupling 4-terminal junction with parameters «', 7' and 8', and phase 
differences fa x and fa 2 on both sides of the junction. Taking fa, fa 2 «; 1, after simple algebra we find 

El2 = - K ' cos f- n' cos f - 87' cos f cos ^ + 86' sin & sin f (B2) 

If «' > 7 > i,5', then the state of the qubits will be dominantly determined by fa and fa. In that case, we can 
ignore <j>\ in all the calculations of the qubit states and the results previously obtained will hold. The value of fa { then 
can be obtained from the equilibrium value of tpi (= fa - n) by 

0>2sin-(^) J ,2sin-(^)«^ (B3 , 
where J t is the current in the SQUID loop of the i-th qubit. The last term in (B2) give a coupling potential 

= (^2) (B4) 
In quantum regime the total Hamiltonian fro the two qubits will therefore be 

n tot = -H 1 +n 2 + n coup (B5) 


are the uncoupled Hamiltonian for the i-th qubit and ri coup is the coupling Hamiltonian respectively. 
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